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On  the  Inclination  of  Satellite  Orbits  About  an 
Oblate  Precesslng  Planet 

Introduction 

It  Is  shown  that  satellite  orbits  about  an  oblate 
precesslng  planet  will  maintain  a  constant  inclination  to 
the  planet's  equator  under  a  certain  condition.  This  condi¬ 
tion  Is  that  the  motions  of  the  satellite's  perlcenter  and 
of  Its  ascending  node  on  the  equator  plane  are  rapid  when 
compared  to  the  precession  of  the  planet.  This  Is  the  case 
for  both  satellites  of  Mars  and  explains  why  It  Is  not  Just 
a  coincidence  that  both  of  these  satellites  have  small  (‘^2*) 
Inclinations  to  the  equator  of  Mars  at  the  present  time. 


1.  It  Is  well-Known  that  In  the  absence  of  satellites, 

an  oblate  spinning  planet  whose  equator  Is  inclined  to  Its 

orbit  plane  will  precess  about  the  normal  to  Its  orbit  plane 

with  period  given  by  Ts2  (P^)*  C  seca?  (1-1).  P  Is  the 

3  P  (C-A) 

period  of  rotation  of  the  planet  and  P'  Is  its  period  of 
revolution  about  the  sun.  C  and  A  are  the  planet's  moments 
of  Inertia  about  an  axis  perpendicular  to  the  equator  plane 
and  an  axis  In  the  equator  plane  respectively.  For  sln^li- 
city  we  assume  axial  symmetry  for  the  planet.  ^  Is  the  angle 
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between  the  planet's  orbit  and  equator  planes,  also  known  as 

the  obliquity.  For  Mars  formula  (1-1)  gives  us  T>2.17>^10^ 

years  where  C-A  is  taken  from  satellite  data.^ 

C 

If  ^  is  the  angle  between  the  planet's  equatorial  and 
orbit  planes  then  the  inclination,  relative  to  the  equator,  of 
a  satellite  orbit  which  is  fixed  in  space  will  vary  as  the 
planet  precesses.  If  is  the  angle  between  the  satellite 
and  planet  orbit  planes  then  the  inclination  would  vary  by  2ff 
and  by  26  toz6  —  ^»  For  Mars  ^>25*’  12'.  However, 
Phobos  and  Oeimos  are  Inclined  to  Mars'  equator  by  less  than 
two  degrees.  Similar  results  hold  for  the  inner  satellites  of 
Jupiter,  Saturn,  Uranus,  and  Neptune.  Since  these  planets  are 
all  precessing  it  is  evident  that  they  must  drag  the  orbit 
planes  of  their  satellites  around  with  them.  Otherwise  the 
low  inclinations  of  these  satellites  would  amount  to  an  un¬ 
believable  coincidence.  Two  possibilities  suggest  themselves 
since  the  equator  of  a  planet  is  defined  both  by  the  spin  axis 
of  the  planet  and  by  the  equatorial  bulge  that  this  spin 
produces.  In  the  first  place,  the  spin  can  affect  the  incli¬ 
nation  through  tides  raised  on  the  planet  by  the  satellite. 
These  tides  transfer  spin  angular  momentum  into  orbital 
angular  momentum  of  the  satellite  or  vice  versa  depending  on 
whether  the  orbital  period  or  the  spin  period  is  longer. 
However,  upper  bounds  can  be  put  on  these  tidal  effects  and 
they  definitely  rule  out  this  possibility  as  an  explanation  of 
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the  inclinations  of  the  satellites.  Ihe  second  possibility  is 
directly  related  to  the  equatorial  bulge  of  the  planet.  It  is 
well  known  that  the  major  effect  of  the  oblateness  of  a  planet 
on  its  satellites'  orbits  is  to  produce  a  secular  motion  of 
the  pericenters  and  of  the  nodes.  In  this  paper  it  is  shown 
that  if  these  motions  are  sufficiently  rapid  then  the  inclina¬ 
tion  of  the  satellite's  orbit  relative  to  the  planet's  equator 
is  unchanged  by  the  planet's  precession.  It  is  then  apparent 
that  in  this  case  if  satellites  are  brought  into  their  planets' 
equator  planes  (either  by  being  formed  there  or  by  some  other 
means)  then  their  inclinations  to  this  plane  will  remain  constant 
as  the  planet  precesses. 


2.  The  notation  which  is  followed  in  this  paper  is  set 
out  below. 

J\  is  the  longitude  of  the  ascending  node. 

1  is  the  orbital  inclination  relative  to  the  planet's 
equator. 

u  is  the  longitude  of  the  perlcenter. 

a  is  the  semi -major  axis. 

e  is  the  eccentricity. 

f  is  the  true  anomaly. 

M  is  the  mean  anomaly.  M=n(t-t,) 

n  is  the  mean  motion,  n^^ where 

T 
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T  la  the  period  of  the  satellite. 

i/  la  the  mean  motion  of  the  planet  about  the  sun. 

la  the  planet's  period  about  the  sun. 
s  Is  the  spin  angular  velocity  of  the  planet. 

R  Is  the  radius  of  the  planet. 

M  Is  the  mass  of  the  planet. 

Explanation  of  Pig. (1-2). 

The  two  coordinate  systems  which  will  be  used  are  drawn 

In  Fig. (1-2).  0  la  their  common  origin. 

The  X,  Y,  Z  system  Is  an  Inertial  system  with  the  X  Y 

plane  coinciding  with  the  planet's  orbital  plane,  i,  A, 

are  the  usual  unit  vectors. 

The  X,  y,  z  set  Is  one  In  which  the  x  y  plane  Is  the 

planet's  equatorial  plane.  Ihe  x  axis  lies  In  the  X  Y  plane 

and  makes  an  angle  of  with  the  X  axis.  Hence, =  & Is  the 

angular  velocity  of  precession  and  jrf  Is  the  planet's  obliquity, 
*  ^  ^ 

1,  j,  k  are  the  unit  coordinate  vectors  of  this  system. 
^=//A=/U(sin$i^  j+cos(rf  ic). 
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Let  r  be  the  radius  vector  to  the  satellite  from  0. 

r=|f|  r»a(l-e* )  (1-2). 

r  1-fecosf 

Let  ^  be  a  unit  vector  in  the  orbit  plane  orthogonal  to 
and  in  the  direction  of  increasing  f. 

Finally,  let 

We  will  need  the  matrix  which  relates  ^  to  1,  3» 

This  is  Just  the  Euler  matrix  and  can  be  written  as  follows > 


\ 

i 

0 

s 

A 

j 

A 

/ 

A 

UJ 

where  the  matrix  A  is  given  below 
'  cosucosA.  -sinusinyvcosi ,  cosusinA-fsinucosncosi, 
-sinucosA  -cosusinrt  cosi , -sinusinA+cosucosACOsi , 
sinisinn  '  -sinicoaa  , 

We  have  set  u=w+f  in  the  above. 


sinusini 

cosusini 

cosi 

(1-3). 


3.  The  equation  for  ^  is  given  in  terms  of  the  perturba- 

3 

tion  force  per  unit  mass  F  as  follows 

dl=  1  /  r  I  CO  s  ( w+f )  l».a]  (1-4). 

dt  na(l-e^)*i  la/ 

In  the  X,  y,  z  frame  If  is  made  up  of  two  parts.  The  first 
is  due  to  the  planet's  oblateness  and  the  second  is  the  Inertial 
force  due  to  the  planet's  precession.  The  oblateness  has 
been  treated  by  many  authors  including  Kozai^  whose  work  we 


-7- 


■hall  refer  to  for  results  at  a  later  stage.  In  the  present 
paper  a  development  of  the  Inertial  force  will  be  carried  out. 

Let  Fx  denote  the  Inertial  force  and  l^o  the  force  due  to 
the  oblateness. 

Hence, 

has  the  familiar  form 

The  dot  denotes  differentiation  with  respect  to  time  In  the 
precesslng  x,  y,  z  coordinate  system.  ]1  Is  the  precession 

^  h  h 

angular  velocity  and  Is  expressed  In  terms  of  the  j# 
unit  vectors  as /l3:yL((sln0f  J+costz^  k).  Finally,  since  we  are 

s 

considering  uniform  precession,  we  can  set^^o  and  we  get 

The  next  step  Is  to  express  rcos(a)'ff  In  terms  of 

trigonometric  functions  of  the  true  anomaly  f.  We  shall  treat 
the  "centrifugal  force"  first. 

but  (?-^)=0.  Hence,  we  have  rco8(W+f )[jfix(^*r)*^]=rco8(U)+f ) 
(ju*r )  .  Using  the  transformation  matrix  (1-3)  we  can 

express  r  and  ^  In  terms  of  1,  J,  k. 

A  simple  algebraic  computation  then  yields t 
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-4  co8(a/+£)(/tx  (yCb‘?)«w]  = 

Wx 

(  2co8^  f co8*w  +28111^ £81  -8ln2£8ln2aj)  8ln2n8lnl 


[( 8ln2£co82a;‘*’C082£8ln2(u)co8^A  sin21 
[( 8ln2£cos2u^co82£8ln2(j) 8ln2^ 

-8in2i!^  r(2cos2fco8^  +28ln2f8in*w-8in2£8in2iJ)8inACoai  ' 
_( 8 1  n2  £co  8  2ci;-f  CO  8  2  £  8l  n2cj)  CO  8ACO  8  2 1 
(1-5) 

Next  we  ahall  treat  the  "Coriolis  force" 

2rcos(w+f  »2rco8(a;+£)[/l* 

2rcos (cj+f )  ‘(rr+rf ^)y ^2rcos (u>*f )  Ju« ( rff-er^  alnff^) 

a<l-e*) 


Where  I  have  made  use  of 

r=a(l-e^)=er^8lnf^ 
r+ecoM  a(l-e* ) 


(1-6) 


Using /i=yu.(sln<zfJ+cos«^k)  we  get  after  an  elementary  computation 
2rco8(w+£)[  = 

'  ( 2cos*f  co8*w  +28in*’  fsin^w  -sin2£8in2w)  sinn]  ( l+ecosf ) 

.  (8in2£co82u>tcos2£cos2  )  J 


sinaf(r3f ) 


(sln2£co8  2(i;4'cos2£cos2u;)8inrL  )  esinf 

(sin2f  sin2(»/-co82fcos2i.;-sln^f  sin*iti)coaACOsi 

co8g((r3f )  r{(sin2fcos2uH-cos2£cos2u>)slnl\  (1-+  ecosf ) 
a(i-e*) 

j (sin2£sin2u;-2cos2fco8*w-2sin^f Bin®u»)3ini}  esinfj 

(1-7). 


4.  At  this  stage  we  will  take  the  time  average  of  ^  due 
^  dt 

to  Pi  over  a  single  orbit  of  the  satellite.  In  order  to  carry 
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out  this  average  we  shall  need  a  formula  relating  f  to  M. 

This  formula  follows  directly  from  the  equations  of  motion. 

1)  r-r6^=-G  Mp  11)  1  d(r*6)»o 

r  r  dt 


11)  Implies  r2de=h=constant.  Since  f  differs  from  e  by  a 
Ht 

constant  angle,  df>h  .  Using  formula  (1-6)  we  have 

,  3t 

er*  slnf f sehslnf  =  i 

STiieiy  irn:^) 

r=ehco8fi*eh*cosf 
afl-eij  a{l-;*Jra 

r6*=h* 


hence  1)  crlves  h^  -eh^cosf  =G  M« 

?3  kTr=5«)r5  -F?" 

h  =a  Mpa(l-e^)  is  the  result  we  are  after. 

df»df  dt=l  M=h  =JS~Jf7aJT^e^ 

Si  3t  3m  n  dt  nr*  n  r* 

dfaa^ jl-e*  using  G  Mp»n*a®.  (1-8), 

dM  ri 


We  shall  denote  time  averages  by  a  bar  over  the  averaged 
c|[uantlty. 


Thus  P«  1  (  FdM. 

■5njo 


(1-9). 


In  performing  the  averages  we  can  change  the  Integration 


variable  to  f  thus  obtaining  P»_ 


/i^j  I?) 


Pdf 


(1-10). 


2JT/r- 

Using  this  formula  for  computing  averages  and  noting  that  r 
Is  a  symmetric  function  of  f  about  f^TT,  we  can  Immediately 
observe  that  all  terms  Involving  an  odd  power  of  slnf  will 
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give  zero  upon  averaging.  This  le  the  only  property  of  these 
averages  that  we  shall  need  aside  from  the  obvious  fact  that 
they  are  functions  of  a  and  e  only.  Dropping  the  terms  which 
average  to  zero  we  can  write t 

5i»  1  ain*«f  /  ( 2r^sinZf8in^  +2r^cos*fcos*a>)  ain2nsini 

dt  na*Tl-e*)tkL7  ( 


2r3f  sin^f  (1+ecosf  )8i  rruj  a  i  nn.+r5fcos2f  (l+ecosf  }co  s  2u) 


er^f  8in2f  sinf  aln2(ucoein.cosii 


5 


er3f sln2f sinf sin2uslni 


(1-11) 
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5.  At  this  stage  we  must  refer  to  the  results  of  the 
paper  by  Kozal.  To  a  first  approximation  the  potential  due 
to  an  oblate  planet  can  be  expressed  as 

V  =  where  A  i*  the  angle  of  lati- 

r  r  3 

tude  measured  from  the  equator.  For  an  axially  symmetric 

(about  a  polar  axis)  planet  A^g3 (C-A)  ®  (1-12) 

2  Mp 

where  C  is  the  moment  of  Inertia  about  the  polar  axis  and  A 
is  the  moment  of  inertia  about  an  equatorial  axis.  In  this 
approximation  the  disturbing  potential 

RgQ  MaAa (l-sin^A  ) .  In  Kozal 's  paper  it  is  shown  that 
of  the  six  orbital  elements  a,  e,  1,(J  M  only  the  last 

three  undergo  secular  (non-periodic)  perturbations.  Further¬ 
more,  the  secular  perturbation  of  M  is  a  trivial  one  and  can 
be  considered  as  defining  a  new  or  'renormalized"  mean  motion. 
The  secular  perturbations  in  UJ  andfx  can  be  expressed  in 
terms  of  n  and  A 2  as  follows t 
W=(*4+  2nAt  ri-5sln^l1t 

^  J 

(1-13). 

JT.  “  -A^nco  si  t 


6.  The  last  step  in  the  proof  Involves  averaging  the 

expression  for  ^  over  a  period  of  time  which  is  long  compared 
dt 

to  the  periods  of  revolution  of  the  perl center  and  node.  It 
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is  here  that  the  assumption  of  a  slow  precession  with  respect 
to  these  periods  la  employed.  This  assumption  allows  us  to 
treat  a,  e  and  1  as  constants  In  this  average  and  only  con¬ 
sider  the  motions  of  uj  and since  the  oblateness  only  produces 
secular  variations  of  these  latter  two  elements  as  mentioned 
In  section  5. 

If  we  exeunlne  the  equation  for  ^  under  the  conditions 

dt 

above,  namely,  that  6U  andj^  are  linear  functions  of  the  time 

and  a,  e  and  1  are  constants,  we  see  that  all  terms  on  the 

right-hand  side  have  zero  time  average  providing  that  we 

average  over  a  sufficiently  long  period  of  time. 

For  the  most  Interesting  case  of  small  Inclination  we 

have  d(ite-2dn.  and  It  is  seen  by  inspection  that  upon  averaging 
dt  dt 

over  a  time  equal  to  T«gr«2Fa^  ( 1-e^  gives  us 

^  nAi 

1  I  dl  dtaO  (1-14).  In  special  cases  such  as  sln*l«f4  we 
T  ^at  5 

have  dM«0  and  In  this  case  terms  which  are  independent  of  Xi 
dt 

must  be  averaged  over  much  longer  periods  of  time  before  they 

give  zero.  Other  pathological  cases  occur  when  da;s-dn.  and  In 

dt  dt 

similar  circumstances.  These  cases  give  rise  to  products  of 
trigonometric  functions  whose  argvunents  change  with  time  at 
approximately  the  same  rates  and  this  also  increases  the  time 
Interval  over  which  these  terms  average  to  zer6~. 
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7.  The  condition  under  which  the  Inclination  le  es¬ 
sentially  constant,  except  for  small  periodic  terms,  Is 
expressed  In  a  more  Illuminating  form  In  this  paragraph. 
From  formula  (1-1)  we  have 


li,=3nj^  (C-A)co8gf 
2s  C 


Al8oA=dA.»  Ain 


when  slnl«l.  Using  AjsS  (C-A)  we  get 

2  Mp 


A=  naCaecg(  a3 /na\  aec^  R 
M  naMpaa  (l-e^iz  fit?*/  (l-e«)  a 


(1-15). 


Nbw  the  condition  for  constant  inclination  becomes 


■^>1 

M 


or  3/ns\  sec^  /R\^»l 


For  Mars  we  have  A-l«25»10®  for  Phobos  while  A^S.Oi^lO^ 

fj.  U 

for  Delmos. 

We  see  from  the  magnitudes  of  for  Phobos  and  Delmos 

tJi- 

that  these  satellites  should  undergo  no  change  In  Inclina¬ 
tion  due  to  the  precession  of  Mars.  At  present  the  orbit 
of  Phobos  Is  inclined- by  l„8*and  that  of  Delmos  by  1.4*  to 
the  equator  plane  of  Mars.  Since  the  obliquity  of  Mars  Is 
25^  12*  It  would  be  a  remarkable  coincidence  If  the  Inclina¬ 
tions  of  these  satellites  varied  by  about  50°  (see  section  1) 
and  Just  at  that  time  we  came  to  observe  them  they  both  were 
Inclined  by  less  than  2°  to  Mars'  equator.  In  light  of  the 
previous  discussion  we  can  see  that  this  Is  not  a  coincidence 
at  all  since  the  Inclinations  of  these  satellites  do  not  vary 
as  Mars  precesses. 
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8.  In  this  Bectlon  a  more  refined  calculation  on  orbits 

of  very  low  eccentricity  e«l  and  small  inclination  slnl«l 

is  carried  out.  The  amplitudes  of  the  periodic  terms  in 

the  inclination  are  evaluated.  Prom  formula  (1-11)  we  have, 

setting  e>o,  and  working  to  first  order  in  sini 

3raii?(-Bln^gfsin2ft.sinl'»eln2f^sinn)*  slni^slnn. 
dt  4n 

Since  4|;«l  for  all  satellites  and  planets  in  the  solar  sys- 
n 

tern  we  can  approximate  2T*>U*in^sinn.  (1-16).  Using  equation 

dt 

(1-13)  we  get 

(c-A)nan--9(c-A)  ^ 
a^  2Ma*  16  C  a* 

Integrating  ^*yiain^slnn.  yields  sin^coarv+ro  where 

r<t  is  the  constant  part  of  T.  Finally,  substituting  for  fi- 
and  A  we  get  r=j'3n^^^|^co8g(jco8n-»r.>  (1-17). 

For  the  satellites  of  Mars  r=3.45>»10"®co8iv+ro  for 
Phobos  and  r=l. 38*10”^  co8j\+ro  for  Deimos  or  r=.71"  cosn.+To 
for  Phobos  and  r=2.84''  cosa+Io  for  Deimos. 


9.  Similar  results  hold  for  Jupiter,  Saturn,  Uranus  and 
Neptune.  In  all  cases  the  major  satellites  of  these  planets 
are  close  enough  to  their  planets  to  maintain  fixed  inclina¬ 
tions  with  respect  to  the  equator  of  the  planet. 

Some  of  the  outer  satellites  of  the  major  planets,  as 
well  as  the  Moon,  have  inclinations  which  remain  constant 


-15- 


with  respect  to  their  planet's  orbit  planes.  In  all  these 
cases,  except  that  of  the  Moon,  the  planets  could  mulntaln 
constant  Inclinations  of  the  satellites  with  respect  to 
their  equators  If  no  forces  other  than  those  of  the  planet 
acted  on  these  satellites.  However,  the  sun  produces  signi¬ 
ficant  perturbations  of  the  elements  of  these  outer  satel¬ 
lites.  In  particular.  If  we  can  neglect  the  oblateness  of 
a  planet  then  the  sun  produces  a  motion  of  the  node  of  a 
satellite.  This  motion  Is  uniform  on  the  planet's  orbit 
plane.  This  Is  the  major  motion  of  the  node  for  some  outer 
satellites  and  makes  Inclinations  constant  with  respect  to 
the  planetary  orbit  planes  even  though  the  planets  are 
oblate  and  precesslng. 

It  should  be  pointed  out  that  the  calculations  for 
these  major  planet-satellite  systems  are  not  quite  as  ele¬ 
mentary  as  those  for  Mars.  This  Is  for  two  reasons.  In 
the  first  place,  mutual  Interactions  between  the  satellites 
produce  additional  motions  of  their  nodes  and  perlcenters 
which  may  be  comparable  to  those  produced  by  the  planets' 
oblateness.  Secondly,  the  torque  on  the  planet-satellite 
system  due  to  the  sun  may  Include  significant  contributions 
from  couples  acting  on  the  satellites  In  addition  to  the 
torque  acting  on  the  planet  itself.  For  this  reason  the 
formula  (1-1)  Is  not  always  applicable. 
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Flnally,  in  the  Earth-Moon  case  It  la  well-hnown  that 
the  Moon  produces  a  greater  torque  on  the  Barth  than  the 
sun  does  so  (1-1)  does  not  apply  here  without  modification 
either. 


On  the  Eccentricity  of  Satellite  Orbits  In  the 
Solar  System 

Introduction 

In  this  paper  the  secular  changes  in  the  eccentricities 
of  satellite  orbits  In  the  solar  system  are  Investigated. 

Two  mechanisms  which  affect  the  eccentricities  are  considered. 
One  of  them  Is  the  tide  raised  on  the  planet  by  the  satellite, 
which  has  been  the  subject  of  discussion  In  the  past;  the 
other  Is  the  tide  raised  on  the  satellite  by  the  planet.  It 
Is  seen  that  cases  arise  In  the  solar  system  In  which  each 
of  these  tide's  effect  on  eccentricity  Is  dominant. 


1.  Darwin^  (1909),  Groves^  (1960)  and  Jeffreys^  (1961) 
have  given  arguments  to  show  that  In  most  cases  the  tide 
raised  by  a  satellite  on  a  planet  tends  to  increase  the 
eccentricity  of  the  satellite's  orbit.  If  we  look  at  the 
values  of  the  eccentricities  which  arise  In  the  orbits  of 
the  Inner  satellites  of  planets,  we  see  that  they  range 
down  as  low  as  10“^  for  Tethys.  Since  it  seems  hard  to 
Imagine  any  process  of  formation  of  satellites  which  could 
produce  Initial  values  of  eccentricity  as  low  or  lower  than 
these.  It  would  appear  necessary  to  look  for  some  mechanism 
that  could  produce  a  secular  decrease  of  eccentricity  which 
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could  rival  in  magnitude  the  secular  change  due  to  the  tides 
raised  on  the  planets.  Such  a  mechanism  was  proposed  by 
Urey^  (1958)  In  the  form  of  tidal  wor)clng  In  the  satellite 
due  to  tides  raised  by  the  planet.  If  we  only  consider  the 
case  where  the  satellite  always  presents  the  same  face  to 
the  planet  (this  Is  the  only  case  which  Is  observed  In  the 
solar  system  and  Includes  the  greater  satellites  of  Jupiter 
and  Saturn  as  well  as  the  Moon.  Jeffreys^  1952)  then  It  Is 
easy  to  see  why  the  tide  raised  on  a  satellite  tends  to 
decrease  Its  eccentricity.  The  eccentricity  e»/l+'2  E  L*  (2-1) 

y  m|m^G 

where  E  Is  the  energy  of  the  orbit,  L  Is  the  angular  momentum 
and  and  Mg  are  the  planet  and  satellite  masses.  If  the 
satellite  Is  not  spinning,  then  the  tide  raised  on  It  can 
only  produce  a  radial  perturbation  force.  This  means  that 
L  Is  not  changed  by  the  tide.  Since  any  energy  dissipation 
In  the  satellite  decreases  E  and  since  we  have  E<0,  0<e<l 
and  L  constant,  we  find  that  e  Is  decreasing  also.  If  e^O 
such  dissipation  must  take  place  since  the  height  of  the  tide 
will  vary  with  the  oscillation  In  distance  between  the  satel¬ 
lite  and  planet. 

In  this  paper  It  will  be  shown  that  Insofar  as  their 
effects  on  eccentricity  are  concerned,  the  tides  on  the 
satellites  are  probably  more  Important  than  tides  on  the 
planets  In  all  cases  where  tidal  effects  might  be  significant 
except  for  Phobos  ,  Delmos  and  possibly  the  Moon  and  Jupiter 
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V.  In  the  cases  of  Phobos,  Delmos  and  Jupiter  V,  however, 
the  tides  raised  in  their  planets  tend  to  secularly  decrease 
their  eccentricities  anyway  so  it  is  not  surprising  to  find 
that  the  orbits  of  these  satellites  have  low  values  of  eccen¬ 
tricity. 


2.  In  this  section  the  following  idealized  case  of  a 

planet  with  a  single  satellite  will  be  considered.  The 
assumptions  are  as  follows* 

i)  The  mass  of  the  satellite  is  neglected  in  comparison  with 
that  of  the  planet. 

ii)  The  inclination  of  the  satellite's  orbit  plane  to  the 
planet's  equator  is  taken  as  zero. 

ill)  The  planet  and  satellite  are  both  considered  to  be 

homogeneous  incompressible  spheres  which  can  be  characterized 

by  two  parameters,//, and  Q.  //,  is  the  rigidity  and  Q  the 

specific  dissipation  function^  Qg2nE*  (2-2) 

lAE  dt 

J  dt 

where  E*  is  the  peak  energy  stored  in  the  system  during  a 
cycle  and^dE  dt  is  the  energy  dissipated  over  a  complete 

J  3E 

cycle.  Q  will  in  general  vary  with  the  frequency  and  eunpll- 
tude  of  the  tide  and  the  size  of  the  sphere  in  addition  to 
its  composition. 

Iv)  Finally,  I  will  work  only  to  first  order  in  eccentricity 
in  the  interest  of  simplicity. 
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The  idealized  problem  that  I  have  set  up  has  bean  solved 

C 

neglecting  tides  raised  on  the  satellite  by  Darwin  (1908) 
and  reviewed  by  Jeffreys  (1961).  The  following  will  be  a 
summary  of  Jeffreys'  paper  with  several  corrections  of  mis¬ 
prints.  Jeffreys'  notation  will  not  be  followed;  instead  we 
shall  use  the  following! 

aasemi -major  axis  of  the  satellite's  orbit, 
eseccentricity  of  the  satellite's  orbit. 

G>gravitatlonal  constant, 
n^satelllte' s  mean  motions^  where 
T«aatellite' s  period.  ^ 

Quantities  pertaining  to  the  planet  will  have  a 
subscript  p  and  those  pertaining  to  the  satellite 
will  have  a  subscript  s. 

Msmass  of  sphere. 

Rsradlus  of  sphere. 

^Bdenslty  of  sphere, 
gssurface  gravity  of  sphere. 

“rigidity  of  sphere. 

Q(l/)“Q  as  a  function  of  frequency  V  . 

U)  “angular  rotation  velocity  of  sphere. 

lsl-*3  are  phase  lags  in  the  periodic  tides 
raised  by  the  tidal  potential. 

These  phase  lags  arise  as  follows t  the  tidal  potential 
acting  is  written  as  a  sum  of  periodic  terms  with  different 
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frequencies.  The  response  of  the  tide  to  any  periodic 
component  of  the  potential  will  be  in  phase  with  the  poten¬ 
tial  only  if  no  energy  is  dissipated  by  the  tide.  If  the 
tide  dissipates  energy,  then  its  phase  will  lag  that  of  the 
potential  (at  least  on  the  average).  2&o/  2S, ,  26^,  and  2Ej 
are  the  lags  of  the  tides  with  frequencies  2u-2n,  2U'-3n, 

2h;-n,  and  3  n  respectively. 

1 

Finally,  using  the  notation  set  up  above,  we  can  state 
Darwin's  and  Jeffreys'  result. 


where  ?  =>d  pM  s,  T  p  denotes  secular  rate  of  change  due 

atat  3it  a[t  _ 

to  tides  on  the  planet  and  d  «  denotes  secular  rate  of  change 

dt 

due  to  tides  on  the  satellite. 

The  formulas  for  A  and  H  in  Jeffreys'  paper,  however, 
contain  several  mistaKes  so  they  will  be  corrected  here. 

The  elastic  tide  raised  by  a  potential  S2  (2-4) 

(where  S2  is  a  spherical  harmonic  of  order  two  and  r  is  the 
distance  from  the  center  of  the  planet)  is  equal  to  ARp  82(2-5) 
at  the  surface  of  the  sphere  where 


^0^  ■‘2 


-  -  -  and  1c9=3  G  Ms  (2-6).  If  we  write  ARn=H.. 

4-air-  ^  P 

then  we  find  that  the  surface  inequality  of  the  planet  pro¬ 
duces  an  external  potential  Uip(r)=3  Q  MpRpHpSt . 

5  r  3 


Using  the  correct  expression  for  H  we  geti 
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.  9(a  Mn)^(3 


top  -  i^ip*|2p*jj.3p]  (s-’: 

|) 


2B?^a9/Up*2g^\)  4-4-2 

If  we  next  conalder  the  tides  raised  on  the  satellite 
we  can  modify  the  above  argument  slightly  to  get  ^ 
must  first  realize  that  tJ^^n.  Hence,  2U^-2nBO, 

2UJ^-3n>B-n,  2(tj|-n>n  and  as  a  consequence  &os>0,  and  £i8*~C2s* 
Furthermore,  Uip  above  Is  the  potential  energy  per  unit 
satellite  mass  due  to  tides  on  the  planet.  Hi  -g^G  Ms  Rg  Hg  S2 

r5 

Is  then  the  potential  energy  per  unit  planet  mass  due  to  tides 
on  the  satellite.  To  change  this  to  potential  energy  per  unit 
satellite  mass,  we  must  multiply  Ui 3  by  ^  .  Using  the  pre¬ 
ceding  results  we  have  from  Jeffreys'  formula  that 

s  .-9  M„(3  Mp/'  (0  M.)i»;p..r „  e. 

3t  •  2  Mga'4t(19yUg+2gg^,R,)  (.IT  ^  J  ^ 


3t  ® 


(2-8) 


>.  3.  Before  we  can  conpare  the  magnitudes  of  ^  and  ^  , 

dt  P  3t 

we  must  relate  the  phase  lags,  2&i,  to  the  specific  dissipa¬ 
tion  function  Q.  This  Is  accomplished  In  the  following 


manner.  Let  W  be  the  work  done  on  one  of  our  homogeneous. 
Incompressible  spheres  by  a  body  force  i{t)  which  Is 
derivable  from  a  potential  UC^).  Then 


-23- 


(2-9)  where  ^("J)  Is  the  velocity  of 
the  material  at  point  r.  Using  the  equation  of  continuity 


we  get  7»(^U)+^  U.  Since  6o  «0.  we 

-  ^  ^  t 

have  dW>  (  yOU^*fid8  (2-10)  where  n  is  the  outward  normal  to 
dt  i^r 

the  surface  of  the  sphere o  If  Uccost/t  and  the  surface 
inequality  is  c  co8(vt-2&)  then  T^*{)x.sin(t^t-2&) .  Hence,  we 
get  dWsKcosvt8in( Vt-2&) .  A  single  integration  gives  to  first 

at 

order  in 


/dW  dtg-Kffsinfig  L  ^ 
at  j/  J  dt 


dt 


E*  = 


I 


dW  dt=K  cos2£.  Using  (2-2)  we  get 

©at  lu 


Q= 


rs*  =  1  or  for  large  Q  2&=1,  (2-11) 

dt  tan2&  Q 


Incidentally,  this  calculation  tells  us  that  S(V)  has  the 
same  sign  as  1/  since 


/dE  dt  <o 

at 


4.  Next  I  will  show  how  Q  compares  for  two  spheres  of 
the  same  material  but  of  different  size.  Ihls  behavior  is 
important  for  explaining  why  the  inner  satellites  of  the 
large  planets  have  circular  orbits,  (e.g.  Saturn  I-*V). 

Qualitatively  one  would  expect  Q  to  increase  with  the 
size  of  the  body  for  the  following  reason.  The  energy 
dissipated,  per  unit  volume,  in  a  cycle,  will  depend  only 
on  the  square  of  the  strains  since  this  frictional  energy 
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dlasipatlon  la  a  local  phenomenon.  The  peak  energy  atored 
per  unit  volume,  on  the  other  hand,  Increaaea  for  a  fixed 
atraln  with  the  alze  of  the  body  alnce  the  atreaaea  are 
Increaaed  over  the  purely  elaatlc  onea  by  the  aelf -gravita¬ 
tion  of  the  aphere.  Hence,  Q  Increaaea  aa  the  bodlea  aelf- 
gravlty  becomea  more  In^rtant  than  Ita  elaatlclty.  The 
relative  Importance  of  elaatlclty  and  aelf-gravlty  of  the 
aphere  entera  Into  the  formula  for  the  tidal  aurface  Ine¬ 
quality  In  the  denominator  (19/i‘f2g^R) .  To  examine  the 
quantitative  dependence  of  Q  on  alze,  we  proceed  aa  followa. 
From  formulas  (2-4),  (2-5),  (2-6)  and  (2-10)  we  have  the 
energy  dlsalpated  In  a  cycle  for  each  periodic  component  of 

the  tide  equals  Ef®  dta  CuR7£>.  where  la 

J  dt  (19«t'*'2cpR) 

Independent  of  fx  and  R.  However,  the  energy  dissipated  must 

be  proportional  to  the  square  of  the  strain  as  mentioned 

previously.  Therefore,  the  energy  dissipated  In  a  cycle 

equals  DlR  where  D|,  Is  Independent  of  u,  and  R  also. 

This  last  formula  Is  derived  as  follows.  Using  (2-5)  and 

(2-6)  we  have  the  surface  strains  proportional  to  R* 

(l^+2g/>R) 

hence  the  square  of  the  strain  Itegrated  over  the  volume  of 

the  sphere  gives  a  result  proportional  to  R^  where 

(i9/t-i>2^R) 

we  assume  the  surface  value  of  strain  throughout  the  sphere 
as  an  approximation  to  get  the  proportional  result.  Com¬ 
paring  these  expressions  we  get  X  where  g  Is 

(19^+2g^Rj 
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Independent  of  ll  and  R.  This  gives  Q  »l4-2anR  (2-12)  where 

Qo  Tf/E" 

Qq  is  the  value  of  Q  for  a  body  where  self-gravity  Is  negli¬ 
gible.  In  the  case  of  liquid  or  gaseous  spheres  and  Qq 
is  not  defined. 


5.  In  this  section  the  two  rates  of  change  of  e  will  be 

compared  both  in  magnitude  and  in  sign. 

First  I  shall  deal  with  the  question  of  sign,  de  „ 

dt 

has  the  sign  of  -(6op-42lE'lp+62B,  ‘♦‘lEao)*  Earth,  Q 

4  4  2  ^ 

and  hence  varies  by  less  than  a  factor  of  four  over  a 

p 

range  of  one  cycle  per  second  to  one  cycle  per  year.  In 

this  case  £i  must  be  the  dominant  term  and  the  sign  of  de  q 

dt  ^ 

is  the  same  as  the  sign  of  2a;-3n.  VThile  this  constant 

behavior  of  Q  with  frequency  may  not  be  true  for  all  planets 

(especially  not  the  major  ones)  it  is  still  likely  that  the 

Si  term  is  dominant  because  of  its  relatively  large  coefficient. 

If  this  &i  term  is  dominant,  we  have  de  for  all  satel- 

dt  ^ 

lltes  except  Phobos,  Deimos,  Jupiter  V  and  the  retrograde 
ones. 


In  the  case  of  tides  raised  on  satellites  the  sign  of 

de  g  is  the  sign  of  -(2^28'*’— *^3s^  where  £2^  and  £3^  have 

dt  2  2  _ 

the  sign  of  n  and  therefore  are  positive.  Hence,  3e  X 0 

dt  ® 

for  all  satellites  which  keep  the  same  face  toward  their 
planets. 
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Next,  we  will  compare  the  magnitudee  of  these  two 
rates  of  change  of  e.  From  formulas  (2-4)  and  (2-8)  we 
get 

s/Rs)  (19//p»2gp/?pRp)|^^2s  ^3aj  (2-13) 

H  p  V^p/  ( i-^Aa+iga/OsRs ijEop”^ip+4^+|iE.3^ 

I  will  examine  this  ratio  In  three  limiting  cases t 

i )  Mp»  2gp/tRp  Zg^gR, 

This  Is  the  case  of  small  satellite  and  small  planet  and 

yields  the  result 

2[e 

3t  *  J25t2,*3e3i  I  (2-14) 

We  see  that  If  the  satellite  Is  appreciably  smaller 
than  the  planet  and  has  approximately  the  same  rigidity  and 
specific  dissipation  function,  we  get  the  tides  raised  on 
the  planet  dominating.  This  case  certainly  applies  to 
Phobos  and  Delmos. 

11 )  /lp«  2gj^pRp  Ms»2g8/>sRs 

This  Is  the  case  of  large  planet  and  small  satellite  and 
we  get 


(2-15) 
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In  this  case  no  general  conclusion  can  be  dra%m  about 
this  ratio. 

ill)  Mp«2gE/>pRp  iL«i« 

'Ihls  is  the  case  of  large  planet  and  large  satellite  and 
gives 


■o 

1 

l^s  /  |£op“i^lp  M 

£+|&3p 

In  this  case  we  see  that  the  satellite  tide  wins.  This 
is  made  even  more  striking  when  we  remember  that  if  the 
ae:elllte  has  the  same  composition  as  the  planet  then 


(2-17) 


6.  A  discussion  of  the  results  obtained  in  the  pro¬ 
ceeding  sections  is  presented.  This  discussion  is  intended 
to  give  an  explanation  of  the  results  obtained  so  far  and 
to  indicate  the  range  of  their  validity. 

The  first  point  dealt  with  will  be  the  effect  that 
tides  raised  on  the  planet  have  on  the  eccentricity.  Let 
us  consider  the  case  where  the  period  of  the  planet's  rota¬ 
tion  is  much  shorter  than  the  period  of  the  satellite's 
revolution  and  the  satellite  is  a  direct  one  (not  retro¬ 
grade).  Since  the  tide  raised  on  the  planet  is  dissipative, 
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\n  have  a  tine  lag  between  the  applied  tidal  force  and  the 
tidal  bulge  It  ralsea.  Because  the  day  is  shorter  than  the 
month  this  time  lag  means  that  the  tidal  bulge  precedes  the 
satellite  in  longitude.  The  effect  of  this  tidal  bulge  lead 
in  longitude  is  to  produce  a  couple  between  the  satellite 
and  planet  which  adds  angular  momentum  to  the  satellite's 
orbit,  at  the  expense  of  the  rotational  angular  momentum  of 
the  planet.  This  is  the  well-)cnown  tidal  couple  which  is 
responsible  for  the  secular  accelaratlon  of  the  Moon. 

So  far  our  argument  has  been  Independent  of  the  ec¬ 
centricity  of  the  satellite's  orbit.  Let  us  next  con¬ 
sider  what  additional  con^llcatlons  arise  when  we  taka  the 
eccentricity  into  account  and  how  they  feed  back  to  effect 
the  eccentricity  which  produced  them.  In  order  to  simplify 
the  picture  let  us  think  of  a  very  eccentric  satellite  orbit. 
The  height  of  the  tide  raised  depends  Inversely  on  the  third 
power  of  the  satellite's  distance  and  the  force  it  produces 
on  the  satellite  Involves  four  more  reciprocal  powers  of 
distance.  Hence,  we  have  the  torque  on  the  satellite  de¬ 
creasing  as  the  sixth  power  of  the  satellite's  distance  from 
the  planet.  This  steep  decrease  with  distance  enables  us  to 
approximate  the  effect  of  the  tidal  bulge  on  the  satellite's 
orbit  by  an  Impulse  at  perlcenter.  With  this  approximation 
the  satellite  must  again  pass  through  the  same  point  at 
perlcenter,  since  bound  orbits  In  inverse  square  force  fields 
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are  periodic.  But,  angular  momentum  has  been  added  to  the 
satellite's  orbit;  hence,  the  apocenter  distance,  and  there- 
force  the  eccentricity,  must  have  been  increased.  In  fact 
the  angular  momentum  of  the  satellite  per  unit  mass  is 
L  »/G  Mpa(l-e*j  =/3  MTpd+e)  where  tp  =  a(l-e)  is  the  dis¬ 
tance  to  pericenter.  If  we  have  AL>0  and  Arp^O  as  discussed 

above,  then  AL^l  go  Ae  is  positive.  Also  from 

2}l  (1+e) 

Arp®0  we  get  A[a(l-e)J  sAa(l-e)-aBeBO  or  Aa^  aAe  so  that 

(l-e) 

Aa  is  positive  also,  as  we  would  expect.  The  previous  dis¬ 
cussion,  when  modified  to  hold  for  smaller  values  of  e, 
accounts  for  the  tendency  of  the  tide  raised  on  the  planet, 
to  increase  the  satellite's  eccentricity. 

The  considerations  presented  above  are  concerned  solely 
with  the  tidal  torque  on  the  satellite.  That  is,  they  only 
make  use  of  the  component  of  the  disturbing  force  which  is 
perpendicular  to  the  satellite's  radius  vector.  In  his 
paper  on  the  Moon's  eccentricity.  Groves  also  considered 
only  the  tidal  torque.  It  is  not  surprising  then,  that  he 
found  the  Moon's  eccentricity  could  only  Increase  due  to  tides 
on  the  Earth.  This  neglect  of  the  radial  components  of  the 
disturbing  force  renders  the  above  arguments,  and  those  of 
Groves  as  well,  Incomplete.  We  shall  consider  how  the  pic¬ 
ture  presented  up  to  now  is  altered  by  the  inclusion  of 
radial  forces. 
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We  ehall  again  take  an  eccentric  orbit  about  our  planet. 
The  relevant  points  are  Identical  with  those  presented  In 
section  1.  There  It  was  shown  that  the  tide  raised  on  the 
satellite  produces  only  radial  perturbation  forces  and  since 
these  cannot  change  the  satellite's  angular  momentum,  but 
must  decrease  Its  energy,  they  must  also  decrease  Its  eccen¬ 
tricity.  TBna  preceding  argximant,  when  applied  to  the  plane¬ 
tary  tide,  shows  us  that  this  tide  may  decrease  as  well  as 
Increase  eccentricity.  The  details  of  whether  we  have  de¬ 
creasing  or  Increasing  eccentricity  depend  on  the  satellite's 
revolution  period,  the  planet's  rotation  period  and  the  euapll- 
tude  and  frequency  dependence  of  Q(t/). 

The  applicability  of  our  results  to  the  actual  planet- 
satellite  systems  extant  Involves  two  questionable  assumptions. 

The  first  assumption  Is  the  neglect  of  all  tides  except 
the  aolld  body  ones.  It  has  recently  been  demonstrated, 
that  In  the  Earth-Moon  system,  the  ocean  tides  which  In  the 
past  were  thought  to  be  of  major  Importance  are  really  much 
less  Important  than  the  solid  body  tides.  This  conclusion 
would  undoubtedly  also  pertain  to  Mars.  For  Jupiter,  Saturn, 
Uranus  and  Neptune,  however,  turbulent  tides  In  their  at¬ 
mospheres  or  possibly  In  any  liquids  which  may  be  found  on 
these  planets,  might  be  of  greater  Importance  than  the  solid 
body  tides.  We  can  still  use  the  two  parameters  yUand  Q  to 
fit  the  tides  on  these  planets,  although  we  can  no  longer 
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hope  to  malce  very  good  estimates  of  their  frequency  and 
amplitude  dependence.  The  satellites  of  these  planets  are 
almost  certainly  solid  since  they  are  not  big  enough  to  have 
held  the  heat  necessary  to  keep  them  liquid  and  are  not  re¬ 
ceiving  enough  heat  to  do  this  either.  Before  leaving  this 
question  of  the  composition  of  the  major  planets.  It  should 
be  mentioned  that  for  Jupiter,  measurements  exist  which  have 
been  used  In  the  past  to  calculate  a  lower  bound  for  Q.  This 
question  will  be  taken  up  In  the  following  paper  where  It  will 
be  shown  that  although  the  measurements  may  be  correct,  their 
Interpretation  is  not. 

The  second  serious  approximation  we  have  made  In  the 

pretense  of  a  linear  superposition  for  the  tides  of  different 

frequencies.  In  developing  formula  (2-7),  Darwin  and  Jeffreys 

both  wrote  the  tide  raising  potential  as  the  sum  of  periodic 

potentials.  They  then  proceeded  to  consider  the  response  of 

the  planet  to  each  of  the  potentials  separately.  At  first 

glance  this  might  seem  proper  since  the  tidal  strains  are 

very  small  and  should  add  linearly.  The  stumbling  block  In 

this  procedure,  however,  is  the  amplitude  dependence  of  the 

specific  dissipation  function.  In  the  case  of  the  Earth,  It 

has  been  shown  by  direct  measurement  that  Q  varies  by  an  order 

of  magnitude  If  we  compare  the  tide  of  frequency  2Ui-2n  with 

the  tides  of  frequencies  2tu-n,  2U)-3n  and  3n.  This  Is  because 

2 

these  latter  tides  have  amplitudes  which  are  smaller  than 
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those  of  the  principal  tide  (of  frequency  2lif-2n)  by  a  factor 
of  eccentricity  or  about  .05.  It  may  still  appear  that  we 
can  allow  for  this  anqplitude  dependence  of  Q  merely  by  adopting 
an  eimplltude  dependence  for  the  phase  lags  of  the  different 
tides.  Unfortunately,  this  is  really  not  sufficient  since 
a  tide  of  small  amplitude  will  have  a  phase  lag  which  Increases 
when  its  peak  is  reinforcing  the  peak  of  the  tide  of  major 
eunplitude.  This  non-linear  behavior  cannot  be  treated  in 
detail  since  very  little  is  known  about  the  response  of  the 
planets  to  tidal  forces,  except  for  the  Earth.  In  our  dis¬ 
cussions  we  shall  use  the  language  of  linear  tidal  theory,  but 
we  must  keep  in  mind  that  our  numbers  are  really  only  para¬ 
metric  fits  to  a  non-linear  problelm. 

There  is  one  more  assumption  which  is  Implied  in  this 
paper.  It  concerns  the  neglect  of  direct  gravitational  in¬ 
teractions  between  bodies  in  influencing  the  eccentricities 
of  the  satellites.  Celestial  mechanicians  in  particular, 

%40uld  consider  this  omission  to  be  a  very  serious  one  over 
periods  greater  than  a  few  thousand  years.  This  is  because 
their  calculations  will  not  guarantee  the  stability  of  sa¬ 
tellite  eccentricities,  perturbed  by  direct  gravitational 
interactions,  for  periods  greater  than  these.  To  this  objec¬ 
tion  one  can  only  offer  the  belief  that  for  well-space  orbits, 
direct  gravitational  interactions  alone  will  not  endanger 
stability  in  eccentricity,  even  over  ages  comparable  to  those 
of  the  solar  system. 
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7.  The  development  of  the  orbit  in  time  will  be  taken  up 

Q 

in  this  section  for  use  later  on.  From  Jeffreys'  paper^  we 
have 

ja  ,  18(G  M^)'^(G  (2-18) 

dt  a'W  (19/ip+2gp^pRp) 

In  the  absence  of  better  information,  Q  is  taken  inde¬ 
pendent  of  frequency  and  amplitude  in  this  section.  This 
enables  us  to  integrate  the  above  equation  which  gives: 

Mp]f^  (G  M)ppRp£p  v/here  a^a  at  t=0. 

L  ll9^p+^gp/)pRpy  J 

If  we  set  &(,  L  ^l  -»3  equal  to  ±  £  then 
5e  ec^G  Mp)  (G  Mb)R^p£)P 

dt  ^  “  a'«19//p+2gp/)pRpy 

de  72  (G  M)  Rs>Os  £>3 

^  a^kd^/ie+ig^Rsl 

oi  is  a  numerical  coefficient  which  depends  on  the  sign  of 
the  various  . 

In  any  case 

5e  _  so/e  \  g/a  \  where  egeo  when  aga^  .  will  be  evalua- 

da  ~  a  \eo/  (^/ 

ted  separately  in  each  case  in  the  following  sections. 


8.  Numerical  estimates  for  ^  are  made  in  this  section 

dt 

for  several  different  satellites. 

1 )  The  Earth-Moon  System 

Since  we  have  more  information  in  this  case  than  in  any 
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Other,  it  will  be  inveetlgated  in  greatest  detail,  Jeffreys^® 
uses  19  >3  for  the  Barth,  with  a  homogeneous  sphere  model. 

If  we  uee  the  seune  value  of  fh  for  the  Moon  as  for  the  Earth, 
then  we  get 


it  *  _  -  28/1738  \/4^  Pop  _  -.7  Pop  where  Q©  i*  the  value  of  Q 

sr:  m^hi  q;; 

IF  P 


for  the  material  when  self -gravity  is  negligible,  as  dis¬ 
cussed  in  section  4.  If  Qos“Qop  then  the  eccentricity  of  the 
Moon's  orbit  would  be  increasing.  However,  there  are  two  con¬ 
siderations  which  tend  to  increase  the  importance  of  the  tides 
on  the  Moon.  In  the  first  place,  the  rigidity  of  the  Moon  is 
likely  to  be  smaller  than  that  for  the  Barth  since  the  high 
rigidity  of  the  Earth  is  due  to  high  pressures  in  the  interior. 
Furthermore,  the  strains  on  the  Moon  are  larger  than  those 
on  the  Earth  by  a  factor  of  4.9,  for  the  value  of  rigidity 
given  above.  It  is  known  that  for  the  Earth  Q  decreases 
with  the  amplitude  of  strain^^  and  this  would  probably  also 
be  the  case  for  the  Moon. 

The  above  would  be  considerably  altered  if  ocean  tides 
were  significant  contributors  to  Q  on  the  Earth,  but  recent 
evidence^^  seems  to  rule  out  this  possibility.  If  we  take 


II 


M.B^then  a  -1.4  Qqp.  We  thus  see  that  small  changes  in 

^  ip 

the  parameters  jj,  and  Q  alter  the  sign  of  ^  so  that  the 

dt> 
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results  for  the  Moon  must  he  consldexed  Inconclusive. 

If  we  determine  the  value  of  1  for  the 

^  )Qp 

Earth,  from  the  present  observed  secular  acceleration  of  the 
Moon,  we  can  use  this  information  to  determine  the  evolution 
of  the  semi-major  axis  of  the  Moon.  MacDonald  has  done  this 
and  claims  that  the  data  is  consistent  with  having  the  Moon 
close  to  the  Earth  between  one-half  and  one  billion  years  ago 
Ignoring  the  difficulties  which  arise  from  such  a  situation, 
we  can  make  the  following  observation  in  regard  to  eccentrl- 
cityt  if  the  rate  of  eccentricity  change  is  dominated  by 
tides  on  the  Earth,  then  we  get  j  ®  (~]^  (2-20) 

For  a^3  10,000  km.,  a  a  380,000  km.  and  e  «  .055  we  get 
e^  a  5*10“®  (Groves  gets  2*10“®)  which  seems  very  small,  es¬ 
pecially  when  compeured  with  the  eccentricities  of  other  sa¬ 
tellites  in  the  solar  systesm,  for  which  planet  tides  produce 
negligible  secular  accelerations.  On  the  other  hand,  if  e 
has  been  decreasing  and  we  assume  its  initial  value  was  less 
than  eo  *  .5,  we  get  the  followings  From  the  result 
^oj  a  j'aj*  we  get  -.6<Jf<0,  so  we  see  that  pj 

and  the  two  rates  of  change  must  have  almost  canceled  each 
other.  In  any  case,  it  seems  likely  from  the  preceding  that 
whether  the  Moon's  eccentricity  is  Increasing  or  decreasing, 
the  two  tidal  effects  are  close  to  being  equal  in  magnitude. 
It  is  worth  noting,  in  this  context,  that  the  Moon's  eccentrl 
city  is  higher  than  that  of  all  other  inner  satellites  in 
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the  solar  system. 

11)  Mars 

As  explained  In  section  5.  Phobos  and  Delmos  are  covered 
by  case  1),  section  5  and  de 


3t  P' 


3e  ^«1. 
dt 

a)  For  Phobos  £2.p<  £>2p'  negative  while 

£3p  Is  positive.  This  tells  us  that  de  p  Is  almost  certainly 

dt 

negative.  This  agrees  with  the  observed  low  eccentricity  of 
.019  for  Phobos.  It  Is  still  necessary  to  show  whether  the 
tides  on  Mars  could  have  appreciably  altered  the  eccentricity 
of  Phobos  In  the  age  of  the  solar  system,  which  we  take  as 
four  billion  years.  Using  the  same  for  Mars  as  Is  used  by 
Jeffreys, we  havet 

3e  -  at  -9.75«10~^5e  sec"^  at  the  present  time.  Next,  using 
at  ^  Qp 

Qa>100,  which  Is  a  typical  value  for  low  eunplltude  tides  on  the 

Earth,  we  get  1  36  ^  -9.75)^10"^^  sec“^.  Since  4x10®  years 

e  at  P 

17  w 

=  1.2X10'^'  sec  we  see  that  e  could  have  been  appreciably  de¬ 
creased  by  tides  on  Mars.  It  should  be  borne  In  mind  that 
the  semi -major  eixls  of  this  satellite  Is  decreasing,  since 
the  satellite's  period  la  shorter  than  the  Martian  day.  This 

means  that  -1  de  was  smaller  In  the  past  than  It  Is  now. 
e  at  P 

Using  formulas  (2-19)  and  (2-20),  one  could  carry  this  analysis 
out  to  Include  the  Integration  over  the  past  four  billion 
years. 

b)  For  Delmos,  we  have  negative  while  Sop,  &2p 

&3p  are  all  positive.  This  assures  us  that  ^  Is  definitely 
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negative  for  Deimos.  Again  using  Jeffreys'  numbers  we  have, 

talcing  Qn  »  100 1  1  3e  „  a  4,5  lO^^O  jgj. -1^  ,q  that  in 

^  i  3t  P 

4x10^  years,  or  1.2  10^^  %ec.,  the  tidal  forces  will  malce  a 
negligible  change  in  e.  For  Deimos,  the  semi-major  eucis  is 
increasing  due  to  the  tides  but  at  such  a  rate  that  this 
effect  may  also  be  safely  neglected. 

Since  the  eccentricity  of  Deimos  is  .003,  we  are  forced 
either  to  accept  the  conclusion  that  Deimos  was  formed 
with  this  initial  eccentricity  or  to  find  some  other  pro¬ 
cess  that  might  have  decreased  the  eccentricity  of  this 
satellite.  In  the  following  paper  one  such  mechanism  is 
proposed.  It  involves  the  direct  gravitational  interaction 
of  Phobos  with  Deimos  coupled  with  the  tidal  forces  on  Phobos. 
It  is  shown  there,  that  the  rate  of  change  of  Deimos'  eccen¬ 
tricity  due  to  this  process,  could  have  been  significant  over 
a  period  of  four  billion  years, 
il)  Jupiter,  Saturn,  Uranus  and  Neptune 

As  discussed  in  section  6,  the  satellites  of  these 
planets  are  likely  to  be  solid,  while  the  state  or  states 
of  the  planets  are  uncertain. 

For  Jupiter  V  the  planetary,  as  well  as  the  satellite, 
tides  decrease  the  satellite's  eccentricity.  These  tidal 
effects  are  very  likely  to  be  significant  and  probably  account 
for  the  satellite's  low  eccentricity  of  .003. 

For  all  the  other  satellites  of  these  planets,  except 
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for  th«  retrograde  ones,  the  planetary  tide  Increases  the 
satellite's  eccentricity.  These  satellite  planet  systems 
probably  can  be  approximated  by  case  11)  of  section  5,  if  we 
assume  /I's  of  the  order  of  those  of  Ice.  If  we  write  equation 
(2-15)  in  terms  of  and  Ogp,  we  gett 


If,  for  simplicity,  we  consider  the  case  of  a  planet  and 
satellite  of  the  same  material,  this  becomes 


Since  ^goPp^D>l .  we  see  that 


all  satellites  with  will  have  decreasing 

eccentricities,  while  those  with  R-<  will  have 

■  iBUg^Rp)* 


Increasing  eccentricities. 

After  this  brief  and  very  speculative  discussion,  we 
can  only  appeal  to  observation,  which  shows  small  eccentrici¬ 
ties  for  the  five  Inner  satellites  of  Jupiter,  the  six  Inner 
satellites  of  Saturn,  the  four  major  satellites  of  Uranus  and 
the  Inner  satellite  of  Neptune.  In  all  cases  where  the  eccen¬ 
tricity  Is  less  than  .01  we  finds  —1  ^  _  ^1. 2)^10^^ sec”^  for 

e  3t 

reasonable  values  of  Q's  and^'s.  This  seems  to  Indicate  that 
tides  raised  on  satellites  are  of  great  significance  In  the 
evolution  of  the  eccentricities  of  these  satellites. 


An  Explanation  of  the  Frequent  Occurrence  of  Near- 
Oommeneurate  Mean  Motions  in  the  Solar  System 

Introduction 

In  this  paper  an  explanation  of  the  improbably  large 
number  of  near-commensurate  pairs  of  satellite  mean  mo¬ 
tions  is  proposed.  It  is  shown  that  special  cases  of 
near-commensurate  mean  motions  are  stable  under  tidal 
forces.  Furthermore,  at  least  four  of  the  best  illustra¬ 
tions  of  commensurabllities  in  the  solar  system  have  this 
stability.  Finally,  the  significance  of  these  stable  con¬ 
figurations  on  the  evolution  of  satellite  systems  is  dis¬ 
cussed. 


1.  The  existence  of  near-commensurabllltes  among  the 
mean  motions  of  the  satellites  and  planets  in  the  solar 
system  has  been  Known  for  many  years.  The  most  famous  of 
these  commensurabllities  Involves  the  Jovian  satellites 
lo,  Europa  and  Ganymede.  Within  observational  accuracy, 
the  mean  motions  (n|  ,  n^  and  n^  respectively)  of  these 
satellites  obey  the  relation  n, -3n2-2r^ -0.  The  motions  of 
these  satellites  have  been  studied  in  great  detail,  first 
by  Laplace,  and  subsequently  by  many  other  authors.  In 
addition  to  this  three-body  commensurablllty,  several  cases 
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of  nttar-coimeniurabilltles  batwaan  tha  maan  motioni  of 
two  aatallitas  hava  alao  baan  known  for  qulta  aoma  tlma. 

Tha  motlona  of  thaaa  palra  of  aatallltaa  hava  alao  baan 
Intanalvaly  atudlad  alnca  thay  ylald  data  from  which  a 
datarmlnatlon  of  tha  aatalllta  maaaaa  can  ba  mada. 

Mora  racantly,  A.  E.  Roy  and  M.  W.  Ovanden^'^  have 
examined  tha  mean  notiona  of  paira  of  planeta  and  aatel- 
litaa  in  a  new  light.  Thay  conaidared  tha  queation  of 
whether  the  obaarved  number  of  near-commanaurata  paira  of 
maan  notiona  in  the  aolar  ayatem  waa  too  great  to  have 
ariaen  from  a  random  diatribution  of  mean  motiona.  Aa 
thia  paper  ia  intended  to  provide  anawera  to  aevaral 
intriguing  queationa  that  they  raiaed,  wa  ahall  begin 
with  a  general  diacuaaion  of  the  contenta  of  their  tvo 
papara. 

In  their  firat  paper,  tha  authora  arrived  at  tha  con- 
cluaion  that  the  preference  for  near-commanaurata  maan 
motiona  in  the  aolar  ayatem  la  inconaiatant  with  the 
aaavimptlon  of  a  random  diatribution  of  mean  motiona  for 
the  planeta  and  aatellitea.  A  aKetch  of  their  proof  of 
thia  important  reault  will  be  preaented  next. 

Before  we  can  prove  anything,  however,  a  aharper  defi¬ 
nition  of  near-commenaurate  mean  motiona  muat  be  given. 

Let  n|  and  n^  (n|>  n^)  be  the  mean  motiona  of  two  bodiea 
about  a  common  center  of  force.  If  two  integer a.  A,  and 

A2  exlat,  auch  that  |n^  -  where  &  la  a  amall 

|ni  All 
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positlve  number,  then  these  mean  motions  are  said  to  be 


nearly  commensurate.  Since  the  ratio  can  always  be  ap- 

ni 

proximated  with  arbitrary  accuracy  by  the  ratio  of  two  in¬ 
tegers,  it  is  necessary  to  limit  the  size  of  the  integers 
considered.  In  their  paper,  Roy  and  Ovenden  arbitrarily  set 
this  limit  for  A|  at  seven.  This  restriction  to  small  in¬ 
tegers  in  no  way  limits  the  scope  of  the  discussion.  In 
fact,  it  can  easily  be  shown  from  perturbation  theory  that 
the  Importance  of  near-commensurabilitles  decreases  as 
their  order  Increases.  Using  our  definition  of  near-commen- 
surability,  we  can  assign  two  integers,  A|  and  A],,  to  every 

pair  of  mean  motions,  n|  and  n^ ,  whose  ratio  n».<l.  Since  the 

ni  T 

smallest  difference  between  adjacent  fractions  is  l-l^l  , 

^  T  H 

there  can  be  at  roost  one  pair  of  integers,  A|  ,A2for  each 


na  -AoJ 
ni  ^ 


«e<l  *.01190. 


pair  of  mean  motions,  nj ,  n2  such  that| 

From  A I  and  A^,  with  A, £7,  we  can  form  17  fractions  with 
valles  from  ^  to  1.  Thus,  given  01190,  the  probability 
that  a  randomly  chosen  ratio  in  the  range  1  to  1  lies  within 

t,of  some  fraction  ^  is  *17>i2e«*  7  »39.67£o  Roy  and 

Ai  ^ 

Ovenden  considered  46  pairs  of  mean  motions  and  compiled  a 
table  which  con^ares  46Pg^,  for  various  01190,  with  the 
observed  number  of  pairs  of  mean  motions  for  which  a  near- 
commensurability  exists  with  &&£«.  This  table,  minus  the 
control  distribution  data,  is  reproduced  below. 
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to 

.0119 

.0089 

.0059 

.0030 

.0015 

46Pe. 

21.7 

16.2 

10.8 

5.3 

2.5 

Observed 
number  of 
pairs  of 
mean  mo- 

33 

26 

20 

12 

6 

tlons  with 
&f  eo 

Table  (3-1) 

Aa  Roy  and  Ovenden  pointed  out,  there  are  two  reasons 
why  the  observed  number  of  near-commenaurabllltles,  which 
are  listed  In  their  table,  might  be  misleading.  In  the 
first  place,  if  n2  has  a  near-commensurabillty  with  both 
n,  and  nj,  then  It  may  also  be  the  case  that  ns  is  nearly 
commensurate  with  n|  .  If  this  is  so,  then  it  is  unclear 
whether  the  commenaurability  between  ni  and  na  should  be 
considered  as  an  independent  one.  In  their  paper,  Roy  and 
Ovenden  showed  that  this  problem  of  "multiple  counts"  was 
likely  to  affect  the  number  of  independent  observed  commen- 
surablllties  listed  in  table  (3-1),  by  2  or  3  for  6*b.0119 
and  even  less  for  smaller  So .  The  second  source  of  error 
arises  from  the  nonuniform  distribution  of  the  ratios 

ni 

on  the  interval  1  to  1.  In  fact,  no  ratio  exists  which 

1 

is  greater  than  .75.  While  it  is  difficult  to  correct 
accurately  for  this  effect,  the  authors  do  show  that  it 
should  not  affect  the  expected  number  of  near-commensurabl- 
lltles,  for  a  given  to  >  ^  ^  factor  very  different  from  .925. 
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In  light  of  the  preceding  discussion,  we  see  that  the 
distribution  of  mean  motions  very  definitely  deviates  from 
randomness  but  that  It  le  difficult  to  say  precisely  how 
large  this  deviation  la. 

In  their  second  paper,  Roy  and  Ovenden  prove  that,  "if 
n  point-masses  are  acted  upon  by  their  mutual  gravitational 
forces  only,  and  at  a  certain  epoch  each  radius  vector  from 
the  (assumed  stationary)  center  of  mass  of  the  system  Is 
perpendicular  to  every  velocity  vector,  then  the  orbit  of 
each  mass  after  that  epoch  Is  a  mirror  Image  of  Its  orbit 
prior  to  that  epoch."  The  authors  call  this  theorem  the 
mirror  theorem  and  the  special  configuration  described  above 
Is  called  a  mirror  configuration.  As  a  corollary  of  the 
mirror  theorem,  the  authors  prove  a  periodicity  theorem 
which  states  that,  "if  n  point-masses  are  moving  under 
their  mutual  gravitational  forces  only,  their  orbits  are 
periodic  if,  at  two  separate  epochs,  a  mirror  configuration 
occurs. " 

After  proving  the  preceding  theorems,  Roy  and  Ovenden 
suggest  that  the  frequent  occurrence  of  mirror  configura¬ 
tions  will  cause  perturbations  on  the  orbits  to  undergo  fre¬ 
quent  reversals  so  that  the  disturbances  they  generate  cannot 
build  up  to  magnitudes  so  large  that  they  endanger  the  sta¬ 
bility  of  the  motion.  For  a  definition  of  stability,  the 
authors  adapt  Poincare's  conditions,  namelyi 
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I)  The  heliocentric  distance  of  any  planet  cannot  increase 
or  decrease  without  limit. 

II)  The  system  repeatedly  passes  through  the  configuration 
It  had  at  time  to,  say  at  times  t|  ,  tj,,  t ^ ,  etc. 

III)  Close  encounters  of  any  pairs  of  planets  are  ruled  out. 
(The  conditions  for  satellite  orbits  are  analogous). 

Finally,  Roy  and  Ovenden  examine  three  of  the  best 

cases  (l.e.,  those  for  which  s  Is  smallest)  of  near-commen- 

surabllltles  In  the  solar  system.  These  Include  three  pairs 

of  satellite  orbits  in  Saturn's  systemi  Hyperion  and  Titan, 

Enceladus  and  Clone,  and  Mimas  and  Tethys.  Ihe  values  of 

na,  -Aa  for  these  satellite  pairs  are  -0.000566,  40.000643 
ni  A| 

and  -0.000784  respectively. 

Observation  provides  the  following  remarkable  results: 
Conjunctions  of  Enceladus  and  Clone  always  occur  near  the 
perlsaturnlum  of  Enceladus.  For  Titan  and  Hyperion,  the 
conjunctions  always  occur  near  the  aposaturnlum  of  Hyperion. 
For  Mimas  and  Tethys,  the  relation  Involves  their  nodes  and 
the  conjunctions  of  these  two  satellles  occur  near  the  mid¬ 
point  between  their  two  ascending  nodes  on  Saturn's  equator 
plane. 

The  examination  of  these  satellite  pairs  thus  reveals 
that  they  satisfy  the  mirror  theorem,  at  least  to  a  first 
approximation.  The  nature  of  this  approximation  and  Its 
significance  will  be  the  subject  of  the  rest  of  this  paper. 
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2.  Before  we  can  proceed  to  the  dlecueslon  of  the  stabi¬ 
lity  of  these  near-commensurate  mean  motions,  we  must  out¬ 
line  classical  perturbation  theory. 

We  shall  describe  the  orbit  of  a  satellite  of  mass  m, 
about  a  planet  of  mass  M,  by  the  following  six  elements t 
a  is  the  semi -major  axis  of  the  orbit, 
e  is  the  eccentricity  of  the  orbit. 

i  is  the  inclination  of  the  orbit  to  the  planet's  equator. 

A  is  the  longitude  of  the  ascending  node. 
t3is  the  longitude  of  the  perihelion. 

tj=/\-hUJ  where  ui  is  the  angle  between  the  ascending 
node  and  the  perihelion. 

&  is  the  mean  longitude  of  the  satellite  at  epoch. 

To  give  the  position  in  the  orbit,  we  will  use  the  mean 
longitude  A  A*  pidt 

n  =/g~R  for  motion  about  a  spherical  planet.  It  is 
a^ 

called  the  mean  motion  since  n*^where  P  is  the  satellite's 

P 

revolution  period. 

For  unperturbed  motion  about  a  spherical  planet,  a, 

e,  i,A ,u7,  and  &  are  constants.  If  the  motion  is  perturbed, 

however,  then  a,  e,  1,a.,D7,  and  e  will,  in  general,  vary 

with  time.  If  VCf)  is  the  total  potential  per  unit  mass 

acting  on  the  satellite,  then  RC^)  is  defined  to  be  vrf)-G  M 

r 

and  is  known  as  the  disturbing  function.  In  terms  of  the 
disturbing  function,  we  can  write  the  equations  of  motion 
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for  a,  9,  and€..  In  what  follows,  R  is  considered 

to  be  a  function  of  a,  e,  l,yx,uT,  and  A.  For  simplicity, 

we  shall  neglect  powers  beyond  the  second  In  the  satellite's 

eccentricity  and  Inclination.  This  allows  us  to  write 

da»2  aR 
3t  na  aX 


deg  -1 

3t  na*e 


dig  -1  aR 
dt  na^l  i/i 


(3-1) 


Perturbations  of  the  first  order  are  obtained  by 

treating  a,  e,  1,  and  n«J(S  M  as  constants  on  the  right-hand 

aAt 

sides  of  the  perturbation  equations,  while  the  mean  longi¬ 
tude,  A  ,  Is  treated  as  a  linear  function  of  the  time.  We 
shall  deviate  from  common  practice,  however,  and  allow  for 
secular  motion  of  GJ  and  when  we  treat  the  periodic  terms 
In  the  disturbing  function.  The  only  restriction  this 
places  on  our  development.  Is  that  It  forces  us  to  treat 
the  secular  part  of  the  disturbing  function  first. 

The  disturbing  function  for  the  action  of  a  satellite 
with  mass  m^  on  one  with  mass  m  Is  given  by 
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R^Gwf  ^  l-X)C^^4y^  422^  j 

Coordinates  are  measured  from  the  center  of  the  planet. 
Primes  refer  to  the  disturbing  satellite. 

Hy-y'9 

It  can  be  shown  that  R  can  be  expanded  in  the  following 
formi^  R"2F(a,  a^ ,  e,  e' ,  1,  1'  )cosT 
where  T«[hA+h^X'  ♦gCT+g'O' ♦fA.+f'jV J  (3-2) 

The  requirement  of  rotational  Invariance  gives  us  the  single 
restriction 

h+h'+g+g' +f+f'»0.  (3-3) 

The  results  obtained  from  first  order  perturbation 

theory  are  only  approximate  due  to  our  treating  a.,  e,  and 

1  as  constants  In  the  right-hand  sides  of  the  perturbation 

equations.  If  necessary,  the  calculations  can  be  extended 

to  higher  order  (the  order  Is  measured  by  the  power  of  nf) 

M 

perturbations.  This  Is  done  by  substituting  the  results  of 
the  first  order  calculations  In  the  right-hand  members  of 
the  perturbation  equations  for  a,  e,  l,/l,u7,  and  g .  We 
shall  need  higher  order  perturbation  theory  when  we  discuss 
commensurablllties  of  more  than  two  satellites.  One  result 
that  we  shall  q[uote  for  later  use  Is  Poisson's  theorem  on 
the  invariability  of  the  semi -major  axis.  This  theorem 
states  that  there  Is  no  secular  term,  due  to  gravitational 
Interactions  between  satellites,  in  the  expression  for  the 
semi-major  eucls,  both  In  the  first  and  second  orders  of 
perturbation  theory. 
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Since  we  will  apply  perturbation  theory  in  cases  where 
near-commensurabllltles  exist,  a  brief  summary  of  the  effects 
of  near-commensurabllities  in  perturbation  theory  will  be 
given  next.  If  we  write  R>2CcosT,  then  the  perturbation 
equations  yield,  upon  integrations 


6;a  ■  2.-  Cl  COST 

(hn+h'  ri  } 

£e  »5  -  CfcosT 

fhn^K^^ 

£1  «  S.-  Ci'cosT 


^0-  CAfinT 

(hn+h'  n'  ) 

cisinT 

Ae=  CasinT 

(hn+h'  n'  ) 


(Note  6|  denotes  first  order  perturbations). 

Here,  in  the  Interest  of  simplicity,  2D  and n  have  been 

treated  as  constants  in  the  right-hand  members  of  the  per- 

t 

turbatlon  equations.  From  the  definition  of  A*  f  ndt-f  £., 
we  see  that 

*  J^^5,ndt+ 

=  C<1'  sinT  +  CV  sinT  ) 

{  (hn+h'  n'  )*  (hn+h'  n'  ) ) 


A  near-commensurabllity  means  that  one  of  the  expressions, 
h*n+l{*  n'  ,  is  very  small  compared  to  n  or  n'  .  From  the 
expression  above,  we  see  that  a  near-commensurabllity  im¬ 
plies  an  enhanced  amplitude  for  perturbations  with  period 
h*n+l{*  o'  .  Since  only  5iA  has  the  small  divisor  squared,  we 
see  that  the  principal  effect  of  a  near-commensurabllity  will 
be  observed  in  the  mean  longitude. 
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3.  We  are  now  in  a  poeltion  to  Investigate  the  stability 
of  near-commensurate  mean  motions.  Our  considerations  will 
find  application  to  the  two-body  cases  of  Enceladus  and 
Dione,  Mimas  and  Tethys,  and  Hyperion  and  Titan,  as  well  as 
to  the  three-body  case  of  lo,  Europa  and  Ganymede.  Other 
possible  exeunples  of  stable  commensurabllitles  to  which  these 
results  might  also  apply  are  mentioned  in  section  10. 

As  a  start,  we  shall  consider  a  planet  surrounded  by 
several  satellites  which  move  on  well-spaced  orbits  of  low 
inclination  and  eccentricity.  We  shall  malce  the  assumption 
that  the  tidal  torques  on  these  satellites  have  produced 
considerable  evolution  of  their  mean  motions  over  a  period 
comparable  to  the  age  of  the  solar  system  (which  we  take  as 
four  billion  years). 

Let  us  make  the  further  assumption  that  the  tidal 
evolution  of  the  mean  motions  of  each  satellite  is  inde¬ 
pendent  of  the  other  satellites.  This  Independent  evolu¬ 
tion  of  mean  motions  is  implied  (at  least  to  second  order 
perturbation  theory)  by  Poisson's  theorem  on  the  invaria¬ 
bility  of  the  semi -major  axis  (see  section  2).  Since,  in 
general,  the  mean  motions  of  different  satellites  will 
evolve  at  different  rates,  the  ratios  of  the  mean  motions 
of  pairs  of  sa tellies  will  vary  with  time.  In  so  doing 
they  will  occasionally  pass  through  a  low  order  commensura- 
billty.  However,  at  any  one  time,  the  ratios  of  the 
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satellltes*  mean  motiona  will  exhibit  a  tendency  for  near- 
commensurability  which  is  consistent  with  a  random  distri¬ 
bution  of  mean  motions.  Such  a  situation  would  certainly 
fall  to  explain  the  strong  tendency  for  near-commensurate 
mean  motions  that  is  observed  among  the  satellites  of  the 
solar  system. 

Zf,  on  the  other  hand,  these  near-commensurabllitles 
were  stable,  then  we  could  account  for  the  large  number  of 
observed  near-commensurablllties.  Suppose,  for  example, 
that  during  the  tidal  evolution  of  mean  motions,  the  ratio 
of  the  mean  motions  of  two  satellites  approaches  very 
close  to  the  ratio  of  two  small  integers.  If  a  near-commen¬ 
surate  motion  of  these  two  satellites  exists,  which  is 
stable  under  further  tidal  evolution,  then  the  satel¬ 
lites  will  remain  in  the  near-commensurabillty  rather  than 
merely  passing  through  it.  However,  the  tidal  torque  on 
each  of  these  satellites  will  not  be  affected  by  such  a 
near-commensurabillty  of  the  satellite  mean  motions.  This 
being  the  case,  in  order  for  the  further  evolution  of  the 
satellite  orbits  to  proceed  without  disrupting  the  near- 
commensurabillty  of  the  mean  motions,  angular  momentum 
must  be  secularly  transferred  between  the  satellites.  At 
first  sight  this  condition  might  appear  to  be  a  violation 
of  Poisson's  theorem  on  the  invariability  of  the  semi -major 
axis.  However,  the  proof  of  this  theorem  Involves  the 
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assumptlon  that  the  arguments  of  the  periodic  terms  in  the 
disturbing  function  are  not  constant.  As  we  shall  see 
later  on,  this  is  Just  the  condition  which  occurs  in  the 
cases  of  stable  commensurabilities. 

V?e  now  see  that  the  orbits  of  a  pair  of  near-coromen- 
surate  satellites  will  still  evolve  as  the  tides  feed 
angular  momentum  from  the  planet's  spin  into  the  satel¬ 
lites'  orbits.  However,  the  satellites  will  share  this 
angular  momentum  between  them  in  just  the  correct  propor¬ 
tion  so  that  their  mean  motions  remain  near-commensurate. 

The  question  of  which  near-commensvurabllitles  are 
stable  will  be  dealt  with  next.  From  the  discussion  of 
the  previous  paragraph  we  see  that  a  necessary  condition  for 
the  stability  of  a  near -commensur ability  is  that  the  direcit 
gravitational  forces  between  the  satellites  involved  are 
strong  enough  to  be  able  to  distribute  the  angular  momentum 
fed  into  the  system  by  the  tides  in  the  manner  necessary  to 
maintain  the  commensurabllity  relation.  Application  of  this 
condition  will  enable  us  to  place  bounds  on  the  tidal  torques 
of  some  satellites.  This,  in  turn,  will  imply  upper  bounds 
for  the  dissipation  of  the  tides  within  the  planets. 

An  examination  of  these  direct  gravitational  forces 
reveals  that  they  decrease  rapidly  as  the  order  of  the  com- 
mensurability  increases  (see  section  9).  This  accounts  for 
the  low  orders  of  the  observed  near-commensxirabilitles  (e.g., 
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2  to  1,  4  to  3  etc.).  When  the  direct  gravitational  forces 
are  so  weak  that  they  cannot  transfer  angular  momentum 
between  the  satellites  at  a  sufficient  rate,  the  two  satel¬ 
lites'  mean  motions  will  evolve  Independently,  each  at  a  rate 
determined  by  the  tidal  torque  on  the  satellite. 

As  discussed  In  section  1,  the  orbits  of  several  pairs 
of  the  best  examples  of  near-commensurabllltles  demonstrate 
remarkable  regularities.  Not  only  do  these  satellite  pairs 
exhibit  near-commensurabllltles,  they  also  show  a  relation 
between  their  conjunctions  and  one  or  more  of  their  orbital 
elements.  Since  the  conditions  relevant  to  the  various  satel¬ 
lite  pairs  differ  only  In  detail,  we  shall  concentrate  our 
attention  on  the  system  of  Enceladus  and  Dlone  whenever  an 
explicit  example  Is  called  for. 

Denoting  the  orbital  elements  of  Dlone  by  primes,  we 
can  state  the  observations  In  the  form  below  (see  section  2 
for  a  definition  of  these  elements). 

2n'  -n-3ro  =  O  and 
dt 

2 a'  -A-<?=V  where  V  oscillates  about 
0^  with  a  small  amplitude.  (Actually,  as  we  shall  see 
later  on,  V  should  oscillate  about  an  angle  very  close  to, 
but  not  equal  to  0” ) .  The  second  relation  states  that  con¬ 
junctions  of  Enceladus  and  Dlone  always  occur  near  the  perl- 
saturnlum  of  Enceladus.  Thus  we  see  that  this  commensurabl- 
11 ty  relation  implies  that  the  term  In  the  disturbing 
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functlon  with  argument  V  can  produce  secular  changes  in  the 
semi -major  axis.  (Actually,  these  secular  changes  only 
occur  when  V  oscillates  about  an  angle  different  from  0°). 
This  Invalidates  the  proof  of  Poisson's  theorem  as  explained 
previously.  It  may  also  be  noted  that  what  we  have  here 
should  not  be  considered  as  a  near-commensurablllty  of  mean 
motions  but  rather  as  an  exact  commensurablllty  Involving 
the  mean  motions  of  Snceladus  and  Dlone  together  with  the 
motion  of  the  perlsaturnlvun  of  Enceladus.  The  fact  that  we 
have  a  near-commensurablllty  of  mean  motions  Is  Just  a  con¬ 
sequence  of  the  small  size  of  ^ in. 

dt/ 

Roy  and  Ovenden  tried  to  show  that  these  near-commen¬ 
surate  satellite  pairs  satisfy  the  hypotheses  of  their 
mirror  theorem.  In  the  approximation  that  the  Inclinations 
of  Enceladus  and  Dlone  are  neglected  and  that  the  eccentri¬ 
city  of  Dlone  Is  taken  as  zero,  we  see  that  this  Is  the  case. 
Furthermore,  even  when  the  eccentricity  of  Dlone  Is  taken 
into  account,  it  may  still  be  argued  that  mirror  configura¬ 
tions  do  occur.  However,  in  the  other  cases  of  near-commen- 
surability  described  by  Roy  and  Ovenden,  mirror  configxira- 
tlons  only  occur  In  a  first  approximation  to  the  actual 
orbits  (e.g.,  when  only  one  eccentricity  Is  taken  as  non¬ 
zero  or  when  both  Inclinations  are  considered  to  be  equal, 
etc.).  Finally,  Roy  and  Ovenden  remarked  that  when  the 
mirror  configuration  was  only  approximately  satisfied. 
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llbratloni  In  the  longitude  of  the  satellites  took  place. 

The  significance  of  these  libratlons  to  their  argument  Is 
unclear  to  the  present  author.  However,  It  Is  well-known 
that  these  libratlons  arise  from  terms  In  the  disturbing 
function  of  very  long  period.  For  Enceladus  and  Dlone,  one 
such  term  has  argument  O'.  This  argument  can  be 

rewritten  as  W=2/\' .  This  Is  Just  the 
term  which  gives  rise  to  the  libra tlon  of  period  3.89  years 
discussed  by  Roy  and  Ovenden  In  their  second  paper. 

We  have  proposed  that  the  near-commensurabllltles  are 
a  consequence  of  the  tidal  evolution  of  satellite  orbits. 

If  this  hypothesis  Is  to  be  tenable,  then  the  tidal  evolu¬ 
tion  of  the  satellite  systems  Involved  must  have  been  ap¬ 
preciable  in  the  age  of  the  solar  system.  Using  the  data 

tabulated  by  Jeffreys^  we  can  express  d£,  ( S  is  defined  by 

dt 

n=noS"^  where  no  is  the  present  mean  motion  of  the  satellite. 

Hence,  dn=-3dS  at  the  present  time)  for  the  various  satel- 

dt  at 

lltes  In  the  solar  system.  In  doing  this,  we  take  a  homo¬ 
geneous  sphere  model  for  the  planet,  fx  Is  Its  rigidity, 

Q  Its  specific  dissipation  function,  yo  its  density,  g  Its 
surface  gravity  and  R  Its  radius.  The  tabulation  below  In¬ 
cludes  all  satellites,  about  planets  having  two  or  more 
satellites,  for  which 


(1+19£^)Q>  10 
2^R 


-17 


-I 


sec 
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Satellite 

/(H-19/1  )Q 
dt  dt  j  2^R 

Phobos 

2.6K10"''' 

Jupiter  V 

T.oao'** 

lo 

4.9X10"'* 

Europe 

1.6x10"''' 

Ganymede 

2.5x10"“ 

Calllsto 

4.0x10"'’ 

Mimas 

4.0x10"''' 

Enceladus 

1.6X10"''' 

Tethys 

3.2X10"'“' 

Dlone 

I.OXIO"'*' 

Rhea 

2.0x10"“ 

Titan 

6.9X10"'* 

Ariel 

2.3x10"“ 

Table  (3-2) 


In  our  notation  1  Is  equivalent  to  sln2&  In  Jeffreys'  nota- 

Q  , 

tlon.  Also  of  Importance  Is  the  fact  that  4>il0^  years 

1.2^10 seconds.  Integrating  Jeffreys'  equation  for 

^  dn,  we  find  that  If  n  Is  the  present  value  of  the  mean 
n  3t 

motion  of  a  satellite,  then  the  mean  motion  of  that  satel¬ 
lite  was  many  times  larger  T  years  ago,  where  T=1  /dl . 

TTlft^ 

^  Is  evaluated  at  the  present  time.  (This  only  applies 

If  ^>0.  If  d|^0,  then  n  was  smaller  In  the  past).  From 
dt  dt 
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thls  formula  we  see  that  for  satellites  whose  present  values 
of  ^  are  greater  than  6.4i^l0"*^  sec”*  ,  appreciable  evolu¬ 
tion  of  their  mean  motions  has  occurred  In  the  age  of  the 
solax  system. 

From  the  discussion  to  follow,  we  shall  see  that  the 
stable  commensur ability  relations  are  between  Mimas  and 
Tethys,  Bnceladus  and  Dlone,  Titan  and  Hyperion,  and  lo, 
Europa  and  Ganymede.  It  Is  also  likely  that  lo  and  Europa, 
Europa  and  Ganymede,  and  Ganymede  and  Callisto  form  stable 
two-body  commensurabllltles  and  that  Dlone,  Rhea  and  Titan 
take  part  In  a  stable  three-body  cominensurablllty. 

From  Table  2  we  see  that  in  all  the  above  mentioned 
commensurabllltles  at  least  one  of  the  satellites  Involved 
has  ^2.69x10"®  sec'*  .  This  indicates  that  even  with 
1  as  small  as  10”**  ,  the  tides  would  still  have 

’Tip  - 

produced  considerable  changes  In  these  systems  In  the  past. 

For  Jupiter  l+19/£_  cannot  be  much  greater  than  1  and  it 
2c^R 

probably  Is  near^l  for  Saturn  also.  This  Implies  that  Q 

|i 

for  Jupiter  and  Saturn  is  about  10^. 

If  we  turn  the  above  argument  around  then  we  can  ask 
the  following  question i  What  percentage  of  all  satellites, 

^  fU  *1 

for  which  dS  /  2 10  sec  and  for  which  another 

at  /  qIi*19U\ 

\  2qoR} 

satellite  exists  with  period  differing  by  less  than  a 
factor  of  seven  from  Its  period,  are  Involved  In  a  commen- 
surablllty  relation?  Inspection  of  Table  2  tells  us  that 
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every  satellite,  except  Phobos,  which  satisfied  the  pre¬ 
ceding  conditions,  is  part  of  at  least  one  definitely 
stable  conmensurabillty  relation.  Furthermore,  it  is  no 
surprise  that  Phobos  is  not  part  of  a  stable  commensurabl- 
11 ty  with  Delmos,  the  other  satellite  of  Mars.  The  mutual 
gravitational  interactions  between  Phobos  and  Delmos  will 
be  very  weak  since  both  satellites  are  extremely  small 
(about  8  and  4  km.  in  radius,  respectively).  If  we  look 
next  at  the  same  group  of  satellites  but  consider  those  for 
which  >  .eaulO"^®  sec*'  then  we  find  that  they 

all  are  involved  in  a  commensurabillty  relation  except  for 
Phobos,  Jupiter  V,  Ariel  and  possibly  Rhea. 

The  arguments  presented  in  this  section  seem  to  imply 
a  tidal  origin  for  the  commensurabllltles  and  a  Q  of  about 
lo'*  for  Jupiter  and  Saturn.  In  the  following  sections  of 
this  paper  other  estimates  of  Q  will  be  obtained. 


4.  In  this  section  the  details  of  a  stability  proof 
for  near-commensurate  satellites  will  be  presented.  As  a 
start,  we  shall  consider  the  following  idealization  of  the 
system  of  Bnceladus  and  Dione.  Two  satellites  with  masses 
m  and  m^  move  in  orbits  about  the  same  planet.  The 
following  restrictions  are  placed  on  their  orbits t 
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in«l  S^«l  *0  i«l' «0  a'>a. 

M  M 


The  quaatlon  of  whether  a  atable  motion  exiati  for  these 

two  hodlea  such  that  2n'-n-^«0,  la  the  one  we  shall  attempt 

Jt 

to  answer.  refers  to  the  observed  secular  rate  of  change 
3t 

of  D. 


In  line  with  the  comment  made  In  section  2,  we  shall 
not  treat  Gas  a  constant,  but  shall  Include  all  motions  of 
VJ  due  to  secular  terms  In  the  dlstvirblng  function  of  m. 
Henceforth,  we  shall  write tus ti|t.  An  Important  distinc¬ 
tion  must  be  made  between  and  ^  Is  the  observed 

2[t  3t 

secular  motion  of  C.  S7|,  however.  Is  Just  that  part  of  this 
secular  motion  which  la  produced  by  secular  terms  In  m's 
disturbing  function.  Due  to  the  commensurabillty  relation, 
the  periodic  terms  In  the  disturbing  function,  which  have 
argument  2/-A-ID',  also  produce  secular  motions  of  ID .  These 
are  Included  In  ^  but  are  not  included  In  tD, . 

In  our  stability  proof,  we  shall  neglect  all  periodic 
temrs  In  the  disturbing  function  except  the  ones  with  argu¬ 
ment  2A'-A-S)’.  It  Is  easy  to  see  why  the  terms  with  this 
argument  might  cause  concern.  By  hypothesis,  this  argument 
has  zero  secular  rate  of  change.  Substituting  these  terms 
In  the  perturbation  equations,  one  sees  that  they  produce 
secular  changes  In  the  elements.  It  might  be  suspected 
that  these  secular  changes  will  disrupt  the  relation 
2n*  -n-S»0.  However,  we  shall  show  that  under  certain 
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clrcunstancea  this  does  not  happen. 

Terms  in  the  disturbing  function  whose  arguments  are 
integer  multiples  of  2A'-A-S^can  be  neglected,  since  their 
coefficients  are  smaller  than  those  of  the  terms  with 
argument  2A'-A-C7by  e  raised  to  the  power  of  the  absolute 
value  of  the  corresponding  Integer  multiple.  All  periodic 
terms  whose  arguments  are  not  multiples  of  2A^  -A -'ll?  produce 
only  small  magnitude,  short  period  perturbations  of  the 
orbits  since  their  amplitudes  are  not  enhanced  by  integra¬ 
tion  (see  section  2).  Finally,  the  theorem  stated  in  sec¬ 
tion  2  concerning  the  invariability  of  the  seml-siajor  axis 
tells  us  that  we  do  not  have  to  worry  about  secular  terms  in 
the  disturbing  function. 

The  stability  proof  runs  as  follows t  Let 
$s2A'-A-C7>  then  the  terms  In  the  disturbing  functions  with 
this  argument  are,  to  first  order  In  e, 

R|  ■-eaB^^(4+fl^)A2(o^)  J  cos  { 
r'i  a-eOm  r(4+#cd  )A2(«*-)1  cos  J 

”5a7‘L  3F  J 

where  ok»a  and  A2(c^}  is  a  Laplace  coefficient.  A2  can  be 
a^ 

expressed  in  terms  of  elliptic  Integrals.  We  shall  need 

the  value  of  C(o^)°(4'»o<d  )A>(o^)  when  we  apply  our  results  to 

doc 

the  system  of  Bnceladus  and  Dlone,  but  at  present  we  shall 
proceed  without  it. 

Next,  the  first  and  second  time  derivatives  of  $  will 


be  evaluated 
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<H-2dA*  -dA-dP»2n '  -n»2d£*  -<ie.-dtS 

3t  dt  3t  at  dt 

||.-2 n'  -n+2 6,'  -  -  S5, »2d£i  -d^ 


In  the  above,  we  have  sett  £*£«‘^&it,  £!•  £,t+ t',t  and 
G7=(7t^(Slt.  £,,  and  CJ|  are  Just  the  motions  of  £,£,' 

and  uj  due  to  the  secular  terms  in  the  disturbing  function. 
Since  we  are  restricting  our  analysis  to  first  order  per¬ 
turbation  theory, 

d£.ade.rsd0r.=0 
dt  dt  dt 

Using  the  perturbation  equations  (3-1),  we  can  now  write t 


dSoane/m/ 
dt 

If-  -n'e^m 


a  >2  dC  (o<)  cos  ^ 
aM  ds*^ 


C(o<.)  +  /a 
a' 


dC 

ast 


cos  $ 


dSLan  ^  m<|^a^|  C(o<.)cos  $ 


(3-4) 


Taking  another  time  derivative  of  f  we  geti 

d*f  sldn' -dn-ne/m'\ /a  dC(oC)sin9dI 

d^  “3t"at  at 

-2n'e/m\  f  C(e<.)+/a  \  dC(oC)l8ingdg 

Im/  L  liT?  dot  J  dt 

»n  lmf\  /a  \  C(Qt)8in6df 
2e\vr}  (av  dt 

Prom  n*a?=G  M,  we  have  dn=-3  /n\da  (3-5) 

dt  Jlaiat 

Using  this  expression  we  can  re-write  the  equation  for  d*| 
as  follows I 
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8lni|S 


(3-6) 


In  some  circumstances,  the  term  containing  sinl 
may  be  neglected.  If  this  term  is  dropped  and  If  we  use  the 
fact  that  C(oi)  Is  positive  (C  will  be  computed  at  a  later 
stage),  then  we  see  that  the  equation  for  $  reduces  to  that 
of  a  simple  pendulum.  As  Is  well-known,  this  equation  can 
always  be  solved  In  terms  of  elliptic  functions.  However, 
for  the  case  of  Bnceladua  and  Dione,  a  small  amplitude 
approximation  is  Justified.  This  yields t 


f  -f  sinyt  where  j/*  =3eC(o<.) 

2 


IS)”"  (1^ 


Hence,  our  result  proves  that  this  special  case  of  a  near- 
commensur ability  will  not  disrupt  itself. 

We  are  now  in  a  position  to  formulate  the  condition 
under  which  the  term  containing  sin  §  ^  may  safely  be 
neglected.  For  simplicity,  we  shall  restrict  ourselves  to 
the  case  where  m^  is  considerably  bigger  than  m.  In  this 
case,  the  ratio  of  the  neglected  term  to  the  one  which  was 
retained  becomes t 


r» 


If  {  «$oSinpt,  then  the  largest  value  that 


r  can  take  is  r*s 


(3-8) 
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The  next  point  to  be  dealt  with  Is  the  stability  of 
the  relation  2n'  -n-^  when  tidal  forces  are  present.  In 
this  connection,  we  will  make  use  of  Jeffreys'  paper  which 
Implies  that  the  rate  of  change  of  mean  motion  due  to  tidal 
forces  Is  given  by 


-dnr=27n^/m\/R 
dt  4  iMlla 


/m\/RVS  2£o 

iMllSrr  (i+ll/u/2g^F 


(3-9) 


where  the  notation  la  that  used  In  section  3.  Rewriting 
equation  (3-6)  to  Include  tidal  effects  (again  we  neglect 
the  sin  S  ^  term),  we  get 


d^  a-  8lnt-»2dnT  -dnr 

d^  “3t  dt 


(3-10) 


If  the  tidal  effects  are  weak  and  If  the  small  angle  approxi¬ 


mation  can  be  used,  this  case  gives 

f  =io8invt+l  /2dnV  -divr\  =  f  8in)^t+  K 
dt  dt  / 


(3-11) 


so  that  the  only  effect  produced  by  the  tides  Is  a  phase 
shift  }  In  $ .  Hence,  we  have  shown  that  weak  tidal  forces 
do  not  upset  the  relation  2n'  -n-^.  The  remarkable  fact 
Is  that  this  relation  Is  maintained  as  the  satellite  system 
evolves  under  the  action  of  the  tides. 


If,  on  the  other  hand,  we  had  a  case  for  which 


2dn(r  -dnr  ,  then  $ig/2dnT  -dnx^  t® 
^[t  l“3t  dt  /  2 


(3-12) 


and  the  relation  2n'  -n-^aQ  would  be  broken  down  by  the  tidal 


forces 
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5.  ‘This  section  will  contain  a 
the  results  of  the  last  section  to 
Bnceladus  and  Olone. 

The  data  Is  given  below. 
Bnceladus 

P-2IE-1. 370218  days 
n 

e  «  .0045 
1  »  0.0« 

3Ss123.43^  per  year 
dt 


numerical  application  of 
the  satellite  system  of 


Dlone 

P'»2n»2.736916  days 
rr 

e'»  .0021 
1'  =  0.0* 

3^^»30.74®  per  year 
dt 


1 


ot  “/a  \  S.631 
W1 


Table  (3-3) 


Using  the  recurrence  relations  for  the  Laplace  coeffi¬ 
cients,  we  arrive  at 

C (o^) «oU 5ok*-2)^b^ ♦( 443o^*-10(X<  )bVa 


where  b^* 


W  I  (l-d^sVe)'^ 

bit  ■  4  f  Afk  -  f^l-c<*8ln*e)*^de] 

T!5[  I  I  J 


Poro<».631  we  get  C(o^)b1.79.  This  establishes  that  C(oC)^0, 
a  fact  which  we  have  used  before. 

Calculations  of  r*  and  can  be  easily  carried  out  now. 


i«*<8.45>iicr^ 


They  yield 
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radians 


(3-13) 


From  the  value  for  r*  we  see  that  our  neglect  of  the 
slni^  term  was  Justified.  Tl'ie  value  of  the  phase  shift 
appears  to  be  quite  small.  Even  if  1 4-19 11=^1 .  it  is  pro- 
bably  too  small  to  be  observed  since  Q  is  likely  to  be 
greater  than  100  for  Saturn  (see  section  3). 


6.  The  system  of  Titan  and  Hyperion  is  similar  to  that 

of  Bnceladus  and  Dione.  If  we  let  primes  denote  quantities 

related  to  Hyperion,  then  we  have  4n'  -3n-(Sf=0  and  the  angle 

dt 

1=4a  -3A-D'  oscillates  about  180“,  This  motion  can  be  shown 
to  be  stable  under  tidal  action  in  exactly  the  same  manner 
as  before.  Unfortunately,  in  this  case,  the  phase  shift 
due  to  tides  is  even  smaller  than  that  for  Bnceladus  and 
Dione  and  offers  no  hope  for  direct  observation.  Finally, 
the  term  in  the  disturbing  functions  with  argument 
4A' -3A-57'  is  Ry=e^Gm'  p7-»o(d  )Aj(o<>)j  cosj^ 


R* 5  ^  'tctd  )  Aj  (ot)J 


cos  $ 


(-14) 


Since  the  coefficient  of  cos  $  can  be  shown  to  be  positive, 
in  this  case,  $  is  stable  about  180^  rather  than  about  0^ 


as  was  the  case  for  Bnceladus  and  Dione 
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7.  In  this  section,  we  shall  treat  the  syetem  of  Mlsiaa 
and  Tethys.  For  this  system,  the  commensurahllity  relation 
Involves  the  node  instead  of  the  perlsaturnlum.  The  con¬ 
junction  of  the  two  satellites  oscillates  about  the  midpoint 
between  their  two  ascending  nodes  on  Saturn's  equator. 

As  an  approximation  to  this  system,  we  shall  neglect 
e  and  e'  and  set  jai^  «  (primes  refer  to  Tethys). 

The  commensurabillty  relation  states  that  §  s4X'  -2A~iV--rL 
oscillates  about  0^.  The  terms  in  the  disturbing  functions 
with  this  argument  can  be  shown  to  have  the  following  format 

R«=-i^an  cos  i 

4a> 

Rjg-i*Gm  (ot)  cos  $ 

4a^ 

where  is  a  Laplace  coefficient  and  is  positive  for 

this  satellite  system. 

Neglecting  the  sin{^  term,  as  in  section  3,  we  get 
d^g-  )/^sin{  where  i^=3i*  Bj  (oc)  ^4n'*  |  m|*n*^m'|  ^a^jj 

For  the  system  of  Mimas  and  Tethys,  a  small  angle 
approximation  cannot  be  used  since  I.  oscillates  about  0* 
with  an  aimplltude  of  95. 3*.  In  this  case,  our  pendultim 
equation  must  be  solved  in  terms  of  elliptic  integrals. 

The  tidal  forces  produce  a  phase  shift  in  §  but  again  we 
can  show  that  it  is  too  small  to  be  observed. 
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8.  The  caae  of  Zo,  Europe  and  Ganymede  was  first  dis¬ 
cussed  by  Laplace  and  since  then  has  been  treated  by  many 
authors.^  In  this  system,  the  commensurablllty  relation 
only  Involves  the  mean  motions  of  these  three  bodies. 

(They  will  be  denoted  by  n,  ,  n^  and  n3).  X|-3Aj.-2A3 
oscillates  about  180^  with  a  very  small  eunplltude.  The 
proof  of  the  stability  of  this  relation,  under  tidal  for¬ 
ces,  Involves  second  order  perturbation  theory;  otherwise, 
it  goes  through  in  exactly  the  seune  manner  as  the  stabili¬ 
ty  proofs  in  the  two  body  cases.  The  stability  for  this 
case  was  known  to  Laplacef  In. this  case  also,  the  phase 
shift  in  1  is  too  small  to  be  observed. 

In  the  past,  many  authors  have  placed  a  lower  limit 3 
on  the  Q  of  Jupiter.  This  was  accomplished  by  noting  that 
observations  of  lo  (Jupiter  I)  gave  no  secular  acceleration, 
to  within  the  observational  accuracy.  Lower  bounds  for  Q 
as  high  as  10*^  have  been  set  by  this  method.  It  is  in¬ 
teresting  to  see  how  the  tidal  stability  of  the  commensu¬ 
rablllty  relation  changes  the  estimates  of  this  lower  bound. 
The  amount  of  angular  momentum  which  the  tides  transfer 
from  the  spin  of  Jupiter  into  orbital  angular  momentum  of 
the  satellites  is  unaffected  by  the  commensurablllty.  How¬ 
ever,  in  order  to  maintain  the  relation  n| -3nx.-2i^  sQ,  this 
anguleu:  momentum  must  be  shared  among  Jupiter's  satellites 
in  a  special  way.  Three  possibilities  will  be  considered. 
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It  is  well-known,  that  in  addition  to  the  three-body 
commenaurablllty  Just  discussed,  Jupiter's  Galilean  satel¬ 
lites  have  the  following  near  two  body  commensurabllitlest 

p 

2nt-n,e0,  2n3-r^  and  7n^-3n3a0.  These  may  or  may  not  be 
stable.  (See  section  9).  If  we  assume  the  stability  of 
the  first  two  two-body  commensurabllltles  and  the  three- 
body  one,  then  the  lower  bound  of  Q  must  be  decreased  by 
a  factor  of  5.  If  we  assume  all  the  above  commensurablll¬ 
tles  are  stable,  then  we  must  decrease  this  bound  by  a 
factor  of  7.4.  If  only  the  three-body  commensurablllty  Is 
stable,  we  again  get  a  reduction  factor  of  about  5. 


9.  The  next  topic  to  be  discussed  will  be  the  evolution 
of  these  stable  commensxurablllty  relations. 

Except  for  the  three-body  case  of  lo,  Buropa  and  Gany¬ 
mede,  the  commensurablllty  relations  we  have  considered  have 
Involved  the  perlcenters  and  nodes  of  the  satellites.  From 
the  stability  equation  (3-10)  for  Enceladus  and  Dlone,  we 
see  that  If  e  Is  below  a  certain  value,  then  the  commensura- 
blllty  relation  Is  broken  down  by  the  tides.  For  Titan  and 
Hyperion,  a  similar  result  holds  for  Hyperion's  eccentricity, 
s'.  The  stability  of  the  Mlmas-Tethys  commensurablllty  would 
be  destroyed  if  1  became  too  small. 


The  Inner  satellites  of  planets  tend  to  have  small 
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eccentrlcititts  and  inclinations.  Since  terms  in  the  disturbing 
functions  with  argument  AiAj  -AjXi  -(A,  -k^)Wi  have  at  least 
A| -A^  powers  of  e|  in  their  coefficients,  and  those  of 
argument  2(A|A2-A2Ai  — (A|-A2)j\()  have  at  least  2(A|-A.2) 
powers  of  1|  in  their  coefficients,  we  see  that  it  is  not 
surprising  that  the  commensurate  relations  we  observe  have 
A|SA2.^1.  In  this  connection,  we  should  note  that,  in  general, 
these  arguments  can  Involve  both  perlcenters  and  nodes,  and 
that  the  only  requirement  on  these  arguments  is  that  of  ro¬ 
tational  invariance,  mentioned  in  equation  (3-3). 

Now  that  we  understand  why  only  commensurabilitles  of 
the  form  A|  where  k  is  a  small  integer,  are  stable 

under  tidal  forces,  we  can  speculate  about  the  past  behavior 
of  satellite  systems.  Appealing  to  the  results  of  the  pre¬ 
vious  paper  we  also  know  that  the  eccentricities  of  many 
satellites  might  have  been  considerably  larger  in  the  past 
than  they  are  at  present.  A  similar  result  has  been  proved 

Q 

for  inclinations  by  Darwin.^  This  suggests  that  in  the 
past,  stable  near-commensurabilitles  might  have  been  much 
more  common  than  at  present,  and  of  considerably  greater 
variety  (i.e.,  larger  differences  A1-A2  were  possible).  Of 
these  commensurabilitles,  all  except  those  which  involved 
the  fewest  powers  of  eccentricity  and  Inclination  would 
then  have  been  disrupted  by  the  tides.  This  would  have 
occurred  as  the  eccentricities  and  inclinations  of  these 
orbits  decreased,  also  due  to  tidal  action.  Remnants  of 


these  higher  order  conmensurabilltles  may  account  for  some 
of  the  near-commensurabllitles  which  we  now  observe,  but  for 
which  no  stability  relation,  such  as  those  discussed,  exists. 

Finally,  there  is  still  another  way  in  which  secular 
changes  in  eccentricity  and  inclination  can  be  produced. 

This  will  be  illustrated  for  the  case  of  Enceladus  and  Dione. 

Referring  to  equations  (3-1),  (3-4),  (3-11)  and  (3-13), 

we  have,  for  the  rate  of  change  of  Enceladus*  eccentricity 

due  to  the  term  in  the  disturbing  function  with  argument  § t 

degQm  C(oc)slnY 
dt  2a 

Since  2r  is  negative,  in  this  case  we  see  that  the  eccentri¬ 
city  of  Enceladus  is  being  decreased  by  the  action  of  Dione. 
This  will  eventually  weaken  the  stability  of  the  commensu- 
rablllty  relation  between  Enceladus  and  Dione  to  a  point 
where  the  tidal  torque  will  disrupt  it.  Possibilities  of 
this  kind  make  it  impossible  to  consider  very  seriously,  the 
exact  present  values  of  the  eccentricities  of  satellites  in 
the  context  of  the  previous  paper  alone.  It  is  also  possible, 
that  a  past  commensurablllty  between  Phobos  and  Delmos  could 
account  for  the  low  value  of  eccentricity  for  Deimos,  as 
discussed  in  the  previous  paper. 

Before  we  end  this  discussion,  it  should  be  mentioned 
that  there  is  a  limit  beyond  which  these  processes  may  not 
decrease  the  eccentricity  of  an  orbit.  That  limit  is  given 
by  the  eccentricity  which  would  be  forced  on  the  orbit  by 
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direct  gravitational  interactions  with  other  satellite*  and 
the  sun. 


10.  In  this  section  further  areas  for  possible  investi¬ 
gation  are  discussed. 

In  the  first  place,  the  stable  commensurabllity  rela¬ 
tions  which  were  discussed  in  this  paper  referred  only  to 
satellite  systems  and  not  to  the  planetary  system.  Tidal 
effects  on  the  planet's  orbits  are  too  small  to  have  any 
significance,  even  over  ages  comparable  to  that  of  the  solar 
system.  However,  the  stability  proof  discussed  in  this  paper, 
would  apply  equally  well  to  other  phenomena  which  might  pro¬ 
duce  secular  changes  in  the  semi-major  axis  of  the  planets. 

In  particular,  during  the  process  of  planet  formation  such 
forces  would  undoubtedly  nave  existed  in  one  form  or  another. 
It  is  then  possible,  that  the  planets  might  also  havfe>  been 
involved  in  commensurabllity  relations  of  the  types  dis¬ 
cussed,  and  that  their  present  distribution  of  mean  motions 
is  at  least  partially  a  reflection  of  these  relations. 

Secondly,  other  stable  commensurahility  relations  may 
exist.  Possible  candidates  appear  to  be  the  two  body  cases 
of  lo  and  Europa,  Europa  and  Ganymede,  and  Ganymede  and 
Calllsto.  For  the  first  two  of  these  pairs,  Griffin^®  has 
remarked  that  the  inner  of  the  pair  is  near  perl Jove  and  the 
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outer  near  apojove,  whenever  a  conjunction  teOcee  place. 
However,  due  to  the  large  perturbations  present  in  these 
systems,  a  much  more  detailed  stability  proof  Is  called  for. 

Finally,  a  possible  three-body  commensurablllty  between 
Titan,  Rhea  and  Olone  is  suggested.  If  we  denote  the  ele¬ 
ments  relevant  to  these  satellites  by  subscripts  3,  2,  and 
1  respectively,  then  the  following  relations  may  be  noted t 
6n5-5n2f2n|  sOOO.081529  degrees  per  day 
Ci +2A^. 08151. 001  degrees  per  day  12. 

Therefore,  6n3-5n2f2n| -lD|-2/^>000.0000±.0009  degrees  per  day 
We  see  that  this  relation  holds  to  within  the  observational 
accuracy  of  six  significant  figures.  If  this  Is  a  stable 
commensurablllty  relation,  then  a  direct  observation  of 
the  tidal  phase  shift  should  be  possible. 
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